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Heavy-ion reactions induced by neutron-rich nuclei provide a unique means to investigate the 
equation of state of isospin-asymmetric nuclear matter, especially the density dependence of the 
nuclear symmetry energy. In particular, recent analyses of the isospin diffusion data in heavy-ion 
reactions have already put a stringent constraint on the nuclear symmetry energy around the nuclear 
matter saturation density. We review this exciting result and discuss its implications on nuclear 
effective interactions and the neutron skin thickness of heavy nuclei. In addition, we also review 
the theoretical progress on probing the high density behaviors of the nuclear symmetry energy in 
heavy-ion reactions induced by high energy radioactive beams. 

PACS numbers: 25.70.-z, 21.30.Fe, 24.10.Lx 
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I. INTRODUCTION 

During the last decade, many radioactive beam facili- 
ties have become available in the world. At these facili- 
ties, nuclear reactions involving nuclei with large neutron 
or proton excess can be studied, providing the opportu- 
nities to study the properties of nuclear matter under the 
extreme condition of large isospin asymmetry. This has 
led to a lot of interests and activities in a new research 
direction in nuclear physics, namely the isospin physics. 
The ultimate goal of studying isospin physics is to extract 
information on the isospin dependence of in-medium nu- 
clear effective interactions as well as the equation of state 
(EOS) of isospin asymmetric nuclear matter, particularly 
its isospin-dependent term, i.e., the density dependence 
of the nuclear symmetry energy. There are already ex- 
tensive reviews on the isospin physics in nuclear physics, 
and they can be found in, e.g., Refs.P, 0, d, 0, [j| @|. 

Knowledge about the nuclear symmetry energy ex- 
tracted from the EOS of isospin asymmetric nuclear mat- 
ter is essential in understanding not only many aspects 
of nuclear physics, such as heavy-ion collisions induced 
by radioactive nuclei and the structure of exotic nuclei, 
but also a number of important issues in astrophysics, 
such as nucleosynthesis during pre-supernova evolution 
of massive stars and the cooling of protoneutron stars. 
Although the nuclear symmetry energy at normal nu- 
clear matter density is known to be around 30 MeV from 
the empirical liquid-drop mass formula @, Q , its values 
at other densities, especially at supra-normal densities, 
are poorly known [l|, Q. Predictions based on various 
many-body theories differ widely at both low and high 
densities 0, [To| . Empirically, the incompressibility of 
asymmetric nuclear matter is essentially undetermined 
even though the incompressibility of symmetric nu- 
clear matter at its saturation density po w 0.16 fm 



has been determined to be 231 ± 5 MeV from nuclear 
giant monopole resonances (GMR) [l| and the EOS at 
densities of 2po < p < 5po has also been constrained by 
measurements of collective flows in nucleus-nucleus colli- 
sions [![• 

Theoretical studies of the EOS of isospin asymmetric 
nuclear matter were started by Brueckner et al. [H[ and 
Siemens [HJ in the late 60's. Since then, There have 
been many studies on this subject based on different 
many-body theories using various two-body and three- 
body forces or interaction Lagrangians. These many- 
body theories provide very useful tools for understand- 
ing the properties of hot and dense nuclear matter, and 
they can be roughly classified into three categories: the 
microscopic many-body approach, the effective-field the- 
ory approach, and the phenomcnological approach. In 
the microscopic many-body approach, the nuclear many- 
body problem is treated microscopically using nucleon- 
nucleon interactions fitted to high-precision experimen- 
tal data and is thus free of parameters. The microscopic 
many-body approach mainly includes the non-relativistic 
Brueckner-Hartree-Fock (BHF) approach [H, [H, [U 13 , 
relativistic Dirac-Brueckner-Hartree-Fock (DBHF) ap- 
proach is na, m i2^j23j23j23 , self-consistent Green's 
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and variational 
In the effective- 
field theory approach, an effective interaction is con- 
structed based on the effective- field theory (EFT), lead- 
ing to a systematic expansion of the EOS in powers 
of density (the Fermi momentum). The effective- field 
theory approach can be based on the density func- 
tional theory 1351 13611 or on chiral perturbation theory 
|42| . Since this approach can be linked 
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to low energy QCD and its symmetry breaking, it has the 
advantage of small number of free parameters and a cor- 
respondingly higher predictive power. The phenomeno- 
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logical approach is based on effective density-dependent 
nuclear forces or effective interaction Lagrangians. In 
these approaches, a number of parameters have to be ad- 
justed to fit the properties of many nuclei. This type of 
models mainly includes the relativistic mean-field (RMF) 
theory [HI El El El ES S3, 51, El JpKrelativis- 
tic and non-relativistic Hartree-Fock [5l|, |52|. l53l l54l l55l . 

I57L 1581 l59l. l6Ct l6ll . l62j or Thomas-Fermi approxima- 
tions [6l|, [sl, |64|] . These phenomenological approaches 
allow the most precise description for the properties of 
finite nuclei. Both the phenomenological and EFT ap- 
proaches contain parameters that are fixed by nuclear 
properties around the saturation density and thus usu- 
ally give excellent descriptions for the nuclear properties 
around or below the saturation density. Their predictions 
at the high density region are, however, probably unre- 
liable. In addition, due to different approximations or 
techniques used in different microscopic many-body ap- 
proaches, their predictions on the properties of nuclear 
matter could be very different even for the same bare 
nucleon-nucleon interaction [H,|6^|. In particular, pre- 
dictions on the properties of isospin asymmetric nuclear 
matter, especially the density dependence of the nuclear 
symmetry energy, are still significantly different for dif- 
ferent many-body theory approaches. 

Fortunately, heavy-ion reactions induced by radioac- 
tive beams provide a unique opportunity to investi- 
gate in terrestrial laboratories the EOS of asymmetric 
nuclear matter, particularly the density dependence of 
the nuclear symmetry energy. During the past decade, 
a large amount of theoretical and experimental efforts 
have been devoted to the study of the properties of 
isospin asymmetric nuclear matter via heavy-ion reac 
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To extract information about the EOS of neutron-rich 
matter, especially the density dependence of the nuclear 
symmetry energy, from heavy-ion reactions induced by 
radioactive beams, one needs reliable theoretical tools. 
Transport models that include explicit isospin-dependent 
degrees of freedom are especially useful for understand- 
ing the role of isospin degree of freedom in the dynamics 
of central nuclear reactions induced by neutron nuclei at 
intermediate and high energies and in extracting infor- 
mation about the EOS of produced neutron-rich mat- 
ter. During past two decades, significant progresses have 
been made in developing semi-classical transport mod- 
els for nuclear reactions. These semi-classical models 
mainly include the following two type s: the Boltzmann- 
Uehling-Ulenbeck (BUU) model [l05j and the quantum 
molecular dynamical (QMD) model jl06l |. While it is 
important to develop practically implementable quan- 
tum transport theories, applications of the semi-classical 
transport models have enabled us to learn a great deal 
of interesting physics from heavy- ion reactions. In par- 
ticular, with the development of the radioactive nuclear 
beam physics, some isospin-dependent transport models 



96], 



m iss IM izn iza iHa isa urn have been successfully 

developed in recent years to describe the nuclear reac- 
tions induced by neutron nuclei at intermediate and high 
energies. 

In studying the properties of asymmetric nuclear mat- 
ter from heavy-ion reactions induced by neutron-rich nu- 
clei, a key task is to identify experimental observables 
that are sensitive to the density dependence of the nu- 
clear symmetry energy, especially at high densities. Be- 
cause of the fact that the symmetry potentials for neu- 
trons and protons have opposite signs and that they are 
generally weaker than the nuclear isoscalar potential at 
same density, most observables proposed so far use differ- 
ences or ratios of isospin multiplets of baryons, mirror nu- 
clei and mesons, such as the neutron/proton ratio of nu- 
cleon emissions [68] , neutron-proton differential flow 
neutron-proton correlation function [8p . £/ 3 He [H 
7I-/7T+ p,|9ill0cLll0l|. E-/E+ [93 and K°/K+ ra- 
tios [l02j |. etc.. In addition, in order to reduce the sys- 
tematical errors, multiple probes taken from several reac- 
tion systems using different isotopes of the same element 
have also been proposed. These multiple probes mainly 
include double ratio or double differential flow. 

Indeed, recent experimental and theoretical analysis 
of the isospin diffusion data from heavy-ion reactions 
has led to significant progress in determin ing the nuclea r 
symmetry energy at subnormal densities [l07l Il08l . Il09l | . 
Based on the same underlying Skyrme interactions as 
the ones constrained by the isospin diffusion data, the 
neutron-skin thickness in 208 Pb calculated within the 
Hartree-Fock ap proa c h is cons istent with available ex- 
perimental data jlldl . Illll 1 1 1 2i | . This symmetry energy 
is also consistent with that from a relativistic mean-field 
model using an accurately calibrated parameter set that 
reproduces both the giant monopole resonance in 90 Zr 
and 20 8 Pb, and the isovector giant dipole resonance of 
208 Pb hlj. It further agrees with the symmetry en- 
ergy recently obtained from isoscaling analyses of iso- 
tope ratios in intermediate-energy heavy ion collisions 



114j . These different studies have provided so far the 



best phenomenological constraints on the symmetry en- 
ergy at sub-normal densities. Information on the symme- 
try energy at supra-normal densities, on the other hand, 
remains inclusive and more efforts are needed to investi- 
gate the supra-normal density behavior of the symmetry 
energy. Heavy-ion collisions induced by future high en- 
ergy radioactive beams to be available at high energy 
radioactive beam facilities will provide a unique oppor- 
tunity for determining the symmetry energy at supra- 
normal densities. 

In the present paper, we review recent progress on the 
determination of the nuclear symmetry energy in heavy- 
ion reactions induced by neutron-rich nuclei. In particu- 
lar, we review the exciting results on density dependence 
of the nuclear symmetry energy at subnormal densities 
determined from recent analysis of the isospin diffusion 
data in heavy-ion reactions. We also discuss the im- 
plications derived from this new information on nuclear 
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effective interactions and the neutron skin thickness of 
heavy nuclei. In addition, we review theoretical progress 
in studying the behavior of nuclear symmetry energy at 
high density from heavy-ion reactions induced by high 
energy radioactive beams. 

The paper is organized as follows. In Section [Til we 
give a brief introduction to the nuclear symmetry en- 
ergy. We then describe in Section HTTI the IBUU04 hadron 
transport model for nuclear reactions induced by radioac- 
tive beams at intermediate energies. In Section IIV1 we 
present the results from the IBUU04 model analysis of 
the isospin diffusion data in heavy-ion reactions and dis- 
cuss the stringent constraint they have imposed on the 
nuclear symmetry energy around the nuclear matter sat- 
uration density. Based on the constrained symmetry en- 
ergy from the isospin diffusion data, we discuss in Sec- 
tion [V] the implications of the isospin diffusion data on 
nuclear effective interactions and the neutron skin thick- 
ness of heavy nuclei. In Section |VT1 we review theoretical 
progress on studying the behavior of the nuclear symme- 
try energy at high density in heavy-ion reactions induced 
by high energy radioactive beams. Finally, a summary is 
given in Section IVlII 



II. THE NUCLEAR SYMMETRY ENERGY 

In the parabolic approximation that has been verified 
by all many-body theories to date, the EOS of isospin 
asymmetric nuclear matter can be written as 



E(p, 5) = E( P , S = 0) + E sym (p)5 2 + 0(5 4 



(1) 



where p = p n + p v is the baryon density with p n and p p 
denoting the neutron and proton densities, respectively; 
S = (Pn—Pp)/(Pp+Pn) is the isospin asymmetry; E(p, 5 = 
0) is the energy per nucleon in symmetric nuclear matter, 
and 
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ld 2 E(p,5), 



8=0 



(2) 



2 8S 2 

is the nuclear symmetry energy. In Eq. ([1]), there are no 
odd-order 5 terms due to the exchange symmetry of the 
proton and neutron in the nuclear matter (the charge 
symmetry of nuclear forces). Higher-order terms in 5 
are negligible. For example, the magnitude of the 5 4 
term at on has been estimated to be less than 1 MeV 
[HI EH, Ell ■ As a good approximation, the density de- 
pendent symmetry energy -E sym (p) can thus be extracted 
from E Bym (p) w E(p, 5=1) — E(p, 5 — 0) which implies 
that the symmetry energy E sym (p) is the energy cost to 
convert all protons in a symmetric nuclear matter to neu- 
trons at the fixed density p. 

Around the nuclear matter saturation density po, the 
nuclear symmetry energy E sym (p) can be further ex- 
panded to second-order as 



^m(p) — E sym (po) + — 



P~ Po 
Po 
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where L and A" sym are the slope and curvature parame- 
ters of the nuclear symmetry energy at po, i- e -, 



L 



3po 



dE sym (p) 
dp 
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K s 



2 d 2 E sym (p) 
ypo g 2p I 



(4) 
(5) 



The L and K sym characterize the density dependence 
of the nuclear symmetry energy around normal nuclear 
matter density, and thus provide important information 
on properties of the nuclear symmetry energy at both 
high and low densities. 

At the nuclear matter saturation density po and around 
(5 = 0, the isobaric incompressibility of asymmet ric nu- 
clear matter can be further expressed as [Ha . Ill5j | 

K{5) = Kq + K aav 5 2 



asy L 



(6) 



(3) 



where Kq is the incompressibility of symmetric nuclear 
matter at the nuclear matter saturation density pp and 
the isospin-dependent part if a sy ~ K sym — 6L [79( char- 
acterizes the density dependence of the nuclear symme- 
try energy. In principle, the information on A" asy can 
be extracted experimentally by measuring the GMR of 
neutron-rich nuclei and a constraint of —566 ± 1350 < 
^asy < 139 ± 1617 MeV has been extracted earlier from 
a systematic study on the GMR of finite nuclei depending 
on the mass region of nuclei and the number of parame- 
ters used in parameterizing the incompressibility of finite 
nuclei [ll[ . The large uncertainties in the extracted -Kasy 
thus does not allow us to distinguish the different nuclear 
symmetry energies in theoretical models. 

Information about the nuclear symmetry energy is im- 
portant not only for nuclear physics, but also for a num- 
ber of important issues in astrophysics. For example, 
the prompt shock invoked to understand the explosion 
mechanis m of a type II supernova requires a relatively 
soft EOS [ll6j , which can be understood in terms of the 
dependence of the nuclear incompre ssibi lity on isospin. 
In the model for prompt explosion [117] . the electron- 
capture reaction drives the star in the latest stage of col- 
lapse to an equilibrium state where the proton concentra- 
tion is about 1/3, which, according to microscopic many- 
body calculations, reduces the nuclear matter incom- 
pressibility by about 30% compared to that for symmet- 
ric nuclear matter. Moreover, the magnitude of proton 
concentration at (3 equilibrium in a neutron star is almost 
entirely determined by the symmetry energy. The pro- 
ton fraction affects not only the stiffness of the EOS but 
also the cooling mechanisms of neutron stars [l!8l Ill9j 
and the possibility of kaon condensation (e - — > K~v e ) 
in dense stellar matter [l20l ] . If the proton concentration 
is larger than a critical value of about 15%, the direct 
URCA process (n — > p + e~ + v e , p + e~ — > n + v e ) 
becomes possible, and would then enhance the emission 
of neutrinos, making it a mor e important process in the 
cooling of a neutron star [118| . 

Unfortunately, the density dependence of E sym (p), es- 
pecially its high density behavior, is poorly known and 
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FIG. 1: (Color online) Density dependence of the symmetry 
energy from the continuous choice Brueckner-Hartree-Fock 
with Reid93 potential (circles), self-consistent Green's func- 
tion theory with Reid93 potential (full line), variational cal- 
culation with Argonne Avl4 potential (dashed line), Dirac- 
Brueckner-Hartree-Fock calculation (triangles), relativistic 
mean-field model (squares), and effective field theory (dash- 
dotted line). Data are taken from [Io| . 



is regarded as the most uncertain among all properties 
of an isospin asymmetric nuclear matter. Even around 
the saturation density, values of the parameters L, K sym , 
and K aay are still very uncertain with different theoret- 
ical models giving very different predictions. This can 
be seen in Fig. [1] where we show the density dependence 
of the nuclear symmetry energy from some of the most 
widely used microscopic many-body theories [Ic| . One 
sees that the theoretical predictions diverge widely at 
both low and high densities. In fact, even the sign of the 
symmetry energy above 3po is uncertain Q. The theoret- 
ical uncertainties are largely due to a lack of knowledge 
about the isospin dependence of nuclear effective inter- 
actions and the limitations in the techniques in solving 
the nuclear many-body problem. 

As mentioned in the Introduction, heavy-ion reactions, 
especially those induced by radioactive beams, provide 
a unique opportunity to pin down the density depen- 
dence of nuclear symmetry energy in terrestrial labora- 
tories. Indeed, significant progress has been made re- 
cently both experimentally and theoretically in determin- 
ing the symmetry energy at subnormal densities. At sub- 
normal densities, a density-dependent symmetry energy 
of E sym (p) » 31.6(/o/po) ' 69 has been found to best repro- 
duce both the isospin diffusion pM 0M [lM EDI, OH 
and isoscaling [114j data in heavy-ion collisions as well as 



the presently acceptable neutron-skin thickness in 208 Pb 
[llOL llllLlll2( ] . Together, these results represent the best 
phenomenological constraints available so far on the sym- 
metry energy at sub-normal densities. Although the high 
density behavior of the symmetry energy remain largely 
undetermined, future high energy radioactive beams to 
be available at high energy radioactive beam facilities 
will allow us to determine the symmetry energy at supra- 
normal densities. 



III. IBUU TRANSPORT MODEL FOR 
NUCLEAR REACTIONS INDUCED BY 
RADIOACTIVE BEAMS 

Transport models are useful theoretical tools not only 
for studying the reaction mechanisms but also for ex- 
tracting information on the properties of produced hot 
dense matter in heavy ion collisions. For nuclear reac- 
tions induced by radioactive beams, comparing experi- 
mental data with transport model calculations allows us 
to extract the information about the EOS of neutron-rich 
matter. The IBUU model, which has been very useful in 
understanding a number of new phenomena associated 
with the isospin degree of freedom in heavy-ion reactions, 
is an isospin- and momentum-dependent transport model 
that is based on the Boltzmann-Uhling-Uhlenbeck equa- 
tion and is applicable for heavy-ion reactions induced by 
both stable and radioactive beams (9(| ■ 

In the IBUU model, besides nucleons, A and N* res- 
onances as well as pions and their isospin-dependent dy- 
namics are included. The initial neutron and proton 
density distributions of projectile and target nuclei are 
obtained from the Relativistic Mean-field model or the 
Skyrme-Hartree-Fock model. It has the option of using 
either the experimental free-space nucleon-nucleon (NN) 
scattering cross sections or the in-medium NN cross sec- 
tions. For NN inelastic collisions, the experimental free- 
space cross sections are used as their in-medium cross 
sections are still very much controversial. The total 
and differential cross sections for all other particles are 
taken either from experimental data or obtained by us- 
ing the detailed-balance formula. Time dependence of 
the isospin-dependent phase-space distribution functions 
of involved particles are solved numerically using the test- 
particle method. In treating NN scattering, the isospin 
dependent Pauli blocking factors for fermions is also in- 
cluded. 

In the following, we outline the two major ingredi- 
ents, i.e., the single- nucleon potential and the NN cross 
sections, of the version IBUU04 of the isospin- and 
momentum-dependent IBUU transport model for nuclear 
reactions induced by radioactive beams [§(|. Other de- 
tails, such as the initialization of the phase space distri- 
butions of colliding nuclei, the Pau li blocking, etc. can 
be found in Refs. 0, 0, M, M, Glzl ■ 
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A. Single-nucleon potential 

One of the most important inputs to all transport mod- 
els is the single-nucleon potential. Both the isovector 
(symmetry potential) and isoscalar parts of this potential 
should be momentum-dependent due to the non-locality 
of strong interactions and the Pauli exchange effects in 
many-fermion systems. In the IBUU04, we use a single- 
nucleon potential derived from the Hartree-Fock approx- 
imation based on a modified Gogny effective interaction 
(MDI) [H3, i-e., 



U(p, S,p,T, x) 



Po Po 

+ B{-?-y (1 - x5 2 ) - Stx— —?—5p Tl 
Po a + 1 pZ 



Po 
2C t ,t 
Po 

Po 



d 3 p' 
d 3 p 



1 + (p - p 1 ) 2 1 'A 2 
1 + (p-ff/A 2 



(7) 



In the above r = 1/2 (—1/2) for neutrons (protons) and 
t / t'; (T = 4/3; f T (r,p) is the phase space distribu- 
tion function at coordinate f and momentum p. The pa- 
rameters A u (x), Ai(x), B, C T;T , C TjT / and A were obtained 
by fitting the momentum dependence of U{p 1 5,p,r,x) 
to that predicted by the Gogny Hartree-Fock and/or 
the Brueckner-Hartree-Fock calculations, the saturation 
properties of symmetric nuclear matter, and the symme- 
try energy of 31 .6 M eV at normal nuclear matter density 
po = 0.16 fm~ 3 |l2l| . The incompressibility Kq of sym- 
metric nuclear matter at po is set to be 211 MeV. The 
parameters A u (x) and Ai{x), given by 



A u {x) 



-95.98- a: - 



2B 



1 



Ai(x) = -120.57 + x- 



2B 



a+1' 

depend on the parameter x that can be adjusted to 
mimic the predicted E syia (p) from microscopic and/or 
phenomenological many-body theories. The last two 
terms in Eq. (|7|) contain the momentum-dependence 
of the single-particle potential. The momentum depen- 
dence of the symmetry potential stems from the different 
interaction strength parameters CV iT < and GV jT for a nu- 
cleon of isospin r interacting, respectively, with unlike 
and like nucleons in the background fields. More specif- 
ically, we use C T:T > = -103.4 MeV and C T:T = -11.7 
MeV. With these parameters, the isoscalar potential es- 
timated from ([/neutron +kproton)/2 agrees reasonably well 
with predictions from the variational ma ny-bo dy theory 
\l22\ , the more advanced BHF approach [l23j | including 
three-body forces, and the Dirac-Brueckner-Hartree-Fock 
(DBHF) calculations [25| in a broad range of density and 
momentum. As an example, we show in Fig. [2] the den- 
sity dependence of the symmetry energy for x = —2, —1, 
and 1. It is seen that the symmetry energy becomes 
softer with increasing value of the parameter x. 

The interaction part of nuclear symmetry energy can 



TABLE I: The parameters F (MeV), G , K syln (MeV), L 
(Me V), an d -Kasy (MeV) for different values of x. Taken from 
Ref. [l09l ]. 



X 


F 


G 


K 


L 


-^asy 


1 


107.232 


1.246 


-270.4 


16.4 


-368.8 





129.981 


1.059 


-88.6 


62.1 


-461.2 


-1 


3.673 


1.569 


94.1 


107.4 


-550.3 


-2 


-38.395 


1.416 


276.3 


153.0 


-641.7 



be parameterized by 

E^(p) = F(x)p/ Po + (18.6 - F(x))(p/p fM 



(9) 



with F(x) and G(x) given in Table U for x = 1, 0, — 1 
and —2. Also shown in Table U are other properties of 
the symmetry energy, including its slope parameter L and 
curvature parameter K sym at po, as well as the isospin- 
dependent part K asy of the isobaric incompressibility of 
asymmetric nuclear matter. It is seen that the stiffness of 
the symmetry energy increases with decreasing x values. 
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FIG. 2: (Color online) Symmetry energy as a function of den- 
sity for the MD I interaction with x = 1, 0, — 1 and —2. Taken 
from Ref. Hoi. 



What is particularly interesting and important for 
nuclear reactions induced by neutron-rich nuclei is the 
isovector (symmetry) potential. The strength of this po- 
tential can be estimated very accurately from (C/ n0 utron — 
c^proton)/2<5 In Fig. El the strength of the sym- 

metry potential for four x parameters is displayed as a 
function of momentum and density. Here we have only 
plotted the symmetry potential at sub-saturation densi- 
ties most relevant to heavy-ion reactions studies at in- 
termediate energies. The momentum dependence of the 
symmetry potential is seen to be the same for all values 
of x as the parameter x appears by construction only 
in the density-dependent part of the single-nucleon po- 
tential given by Eq.®- Systematic analysis of a large 
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MDI with x=1 MDI with x=0 




FIG. 3: (Color online) Sym metry potential as a function of momentum and density for MDI interactions with x = 1, 0, —1 and 
-2. Taken from Ref. [jo^ . 



number of nucleon-nucleus scattering experiments and 
(p.n) charge exchange reactions at beam energies below 
about 100 MeV has shown that the data can be well de- 
scribed by the parametrization U^ lFl „p. = a — bE^ with 
a ~ 22 - 34 MeV and b ~ 0.1 - 0.2 pi fl25l . iHft Imj . 
Although the uncertainties in both parameters a and b 
are large, the symmetry potential at po, i.e., the Lane 
potential, clearly decreases approximately linearly with 
increasing beam energy E^ ul . This provides a stringent 
constraint on the symmetry potential. The potential in 
Eq. ([7]) at po satisfies this requirement very well as seen 
in Fig. [3] This can be more clearly seen from the solid 
line in Fig. 3] which gives the kinetic energy dependence 
of ([/ n cutron — ^proton )/2<5 at normal nuclear saturation 
density given by the MDI interaction with x — 0. Also 
shown in this figure are the predicted kinetic energy de- 
pendence of the symmetry potential from the MDI in- 
teraction with x = for densities away from normal nu- 
clear density, which are presently not known empirically. 
Experimental determination of both the density and mo- 
mentum dependence of the symmetry potential is thus of 
great interest, and heavy- ion reactions with radioactive 
beams provides a unique tool to explore this information 
in terrestrial laboratories. 



Although the effective mass of a nucleon in nuclear 
matter depends the density of nu clear mat ter as well as 
the momentum of the nucleon momen- 
tum dependence of the symmetry potential further leads 
to different effective masses for neutrons and protons in 
isospin asymmetric nuclear matter, i.e., 

rK = ( + rn L dU L \ 
m T I p dp J 

When the effective mass is evaluated at the Fermi mo- 
mentum p T = pf(t), Eq. (flQ|) yields the Landau mass 
which is r elated to the /i Landau parameter of a Fermi 
liquid [23lll28l . ri29l ]. A detailed discussion ab out d ifferent 
kinds of effective masses can be found in Ref. |128| . With 
the potential in Eq. ((7|), the nucleon effective masses 
are independent of the x parameter as the momentum- 
dependent part of the nuclear potential is independent of 
the parameter x. 

Shown in Fig. [5] are the effective masses of neutrons 
and protons at their respective Fermi surfaces as func- 
tions of density (upper window) and isospin asymme- 
try (lower window) obtained from the MDI interaction. 
It is seen that the effective mass of neutrons is larger 
than that of protons and the splitting between them in- 
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FIG. 4: (Color online) Kinetic energy dependence of 
(^neutron — ^proton ) /2S at different densities using the MDI 
interaction with x = 0. The shaded region indicates the ex- 
perimental constraint. 



creases with both the density and isospin asymmetry of 
the medium. Although the momentum dependence of the 
symmetry potential and the associated splitting between 
the neutron and proton effective masses are still highly 
controversial among different approaches and /or u sing 
different nuclear effective interactions [H, Il30l . Il3l| , the 
results presented here are consistent with the predic- 
tions fro m all non-r elativistic microscopic models, see, 
e -g-; [13; Il23t Il32j |. and the non-relativistic limit of 
microscopic relativistic many-body theories, see, e.g., 
[IH HH, [III- Recent transport model studies have in- 
dicated that the neutron/proton ratio at high transverse 
momenta and/or rapidities is a potentially useful probe 
of the splitting between the neutron and proton effective 
masses in neutron- rich matter |90] | . 

The effect due to the momentum dependence of nu- 
clear mean-field potential can be studied by comparing 
its predictions with those obtained using the momentum- 
independent nucleon potential U(p,S,r) = Uo(p) + 

fsym 1 \Pi S, t) with the isoscalar part Uq(p) taken from 
the original momentum-independent soft nuclear poten- 
tial with K = 200 MeV (SBKD) introduced by Bertsch, 
Kruse and Das Gupta [l33j, i.e., 

Usbkd(p) = -356p/p + 303(p/p ) 7 / 6 . (11) 

For the momentum-independent symmetry po- 
tential JJ^ro (Pj S, r), it can be obtained from 



rr MDI(x), c \ 

U sym {p, 6, t) = 
dependent part of 



'sym/dp T using the isospin- 
thc potential energy density 

-symvfv P ■ <> 2 where E syL(p) & given by 
Eq. from the MDI interaction. Therefore, the 

momentum-inde pend ent SBKD potential that has 
Kq = 200 MeV [l33| | and exactly the same E syni (p) as 
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FIG. 5: (Color online) Neutron and proton effective masses 
in asymmetric matter as a function of density (upper win- 
dow) and isospin asymmetry (lower window). Taken from 
Ref. [i07j . 



the MDI interaction is 
^SBkd(p, 5, t) = 



-356 p/p + 303 (p/p ) 7/6 
4rE^(p) + (18.6 -F(x)) 



x {G(x)-l)(p/p f^S 2 



(12) 



B. In-medium nucleon-nucleon cross sections 

Another important quantity in the IBUU model is the 
in-medium NN cross sections. While much attention 
has been given to finding experimental observables that 
can constrain the symmetry energy, little effort has been 
made so far to study the NN cross sections in isospin 
asymmetric nuclear matter. Most of existing works on 
in-medium NN cross sections have concentrated on their 
density d ependence in i sospin symmetric nuclear matter , 

see, e.g., [MliMIMImlmliilliSliMliilliS- 

One simple model for NN in-m edium cross se ctions is 
the effective mass scaling model [134 1 1351 1 1 38j j . In this 
model, while both the incoming current in the initial 
state and the level density of the final state in an NN 
scattering depend on the effective masses of colliding nu- 
cleons, the scattering matrix elements are assumed to 
be the same in free-space and in the medium. As a re- 
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suit, the ratio of the NN cross section in nuclear medium 
! 1 , ! " n to that in free-space value <r N ° N is simply given 



a 



NN 



by 



R 



medium 



^medium 



'NN 



1° 



free 
NN 



(13) 



where /inn and fi* NN are the reduced masses of scat- 
tering nucleon pairs in free-space and in the medium, 
respectively. Since the nucleon mass becomes smaller 
in nuclear medium, the in-medium NN cross section is 
smaller than its value in free space. The relation given 
in Eq. (fl3|) was recently found to be co nsist ent with cal- 
culations based on the DBHF theory 144j for nucleon 
pairs with relative momenta less than about 240 MeV/c 
and in nuclear matter with densities less than about 2pn. 
This finding thus lends a strong support to the effective 
mass scaling model of in-medium NN cross sections in 
this limited density and momentum ranges. We have 
thus extended this model to determine the in-medium 
NN cross sections in the asymmetric nuclear matter pro- 
duced in nuclear reactions induced by radioactive beams. 
For nucleon-nucleon scatterings at higher energies, inelas- 
tic reaction channels become important. Although there 
were some stu dies on in-medium effects in these channels 
|l45l Il46l Il47j | , the experimental free-space cross sections 
are used in the IBUU model as the model is mainly for 
heavy ion collisions at intermediate energies where NN 
inelastic scatterings are less important than elastic scat- 
terings. 

While the effective masses and the in-medium NN cross 
sections have to be calculated dynamically in the evolving 
environment created during heavy-ion reactions, it is in- 
structive to examine the in-medium NN cross sections in 
isospin asymmetric nuclear matter at zero temperature. 
In this situation, the integrals in Eq. ([7]) can be analyti- 
cally carried out. Specifically, it is given by jl2ll . Il48| . 



d 3 p' 



P}(r) 



A 2 



2pA 



■r_ ln (p + Pf(r)) 2 +A 2 



(p-pf(r)) 2 + A 2 



+ 



2p/(r) 
A 



2tan -iP±^_ 2tairl ^M 



A 



A 



(14) 



The reduction factor i? me dium for in-medium NN cross 
sections can thus also be obtained analytically, albeit 
lengthy. 

As an illustration, we show in Fig. [6] the reduction fac- 
tor i? me dium for a simplified case of two colliding nucle- 
ons having the same momentum p. The i? m edium factor 
is examined as a function of density (top panel), isospin 
asymmetry (middle panel) and the momentum (bottom 
panel) . It is interesting to note that not only in-medium 
NN cross sections are reduced compared with their free- 
space values but also the nn and pp cross sections be- 
come different while their free-space cross sections are the 
same. Moreover, the difference between the nn and pp 
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FIG. 6: (Color online) The reduction factor of the in-medium 
nucleon-nucleon cross sections compared to their free-space 
values as a function of density (top panel) , isospin asymmetry 
(mid dle p anel) and momentum (bottom panel). Taken from 
Ref. [l03l. 



scattering cross sections grows in more asymmetric mat- 
ter. The larger in- medium cross sections for nn than for 
pp are completely due to the positive neutron-proton ef- 
fective mass splitting with the effective interaction used. 
This feature may serve as a probe of the neutron-proton 
effective mass splitting in neutron-rich matter. We also 
note that in-medium NN cross sections are also indepen- 
dent of the parameter x as they are solely determined by 
the momentum dependence of the nuclear potential used 
in the model. 



IV. CONSTRAINING THE SYMMETRY 
ENERGY AT SUB-NORMAL DENSITIES USING 
ISOSPIN DIFFUSION DATA 

The IBUU model outlined above have been used to 
study many observables in heavy-ion reactions induced 
by both stable and radioactive beams. In this section, we 
illustrate one of its applications in extracting the density 
dependence of nuclear symmetry energy from the isospin 
diffusion data in heavy-ion collisions. 

Isospin diffusion in heavy-ion collisions has been shown 
to depend sensitively on the densi ty d ependence of 
nuclear symmetry energy [9l], [53, Il49j |. Within a 
momentum-independent transport model, in which the 
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nuclear potential depends only on local nuclear density, 
the isospin diffusion data from recent experiments at the 
NSCL/MSU (National Superconducting Cyclotron Lab- 
oratory at Michigan State University) was found to fa- 
vor a quadratic density dependenc e for the interaction 
part of nuclear symmetry energy [l08| . This conclu- 
sion has stimulated much interest because of its impli- 
cations to nuclear many-body theories and nuclear as- 
trophysics. However, the nuclear potential acting on a 
nucleon is known to depend also on its momentum. For 
nuclear isoscalar potential, its momentum dependence is 
well-known and is important in extracting information 
ab out the equation of state of symmetric nuclear matter 
|, EM, EM, EM, EM, EM, EM, EH [H3 . The momen- 
tu m dependence of the isovector (symmetry) potential 
@i Il2ll . Il26l . Il57l | has also been shown to be important 
for understanding a number of isospin related phenomena 
in heavy-ion reactions [8^, [9(1 HH • It is therefore neces- 
sary to include momentum dependence in both isoscalar 
and isovector potentials for studying the effect of nuclear 
symmetry energy on isospin diffusion. 

Isospin diffusion in heavy ion collisions can in principle 
be studied by examining the average isospin asymmetry 
of the projectile-like residue in the final state. Since re- 
actions at intermediate energies are complicated by pree- 
quilibrium particle emission and production of neutron- 
rich fragments at mid-rapidity, differences of isospin dif- 
fusions in mixed and sym metr ic systems are usually used 
to minimize these effects [108| | . To study isospin diffusion 
in 124 Sn + 112 Sn reactions at E = 50 MeV/nucleon and 
an impact parameter of b = 6 fm, we thus also consider 
the reaction systems 124 Sn + 124 S and 112 Sn + 112 S n and 
at same energy and impact parameter as in Ref. |l08j . 
The degree of isospin diffu sion in the reaction 124 Sn + 
112 Sn is then measured by [1581 ] 

2Xl24g n+ 112 Sn - Xl24 Sn+ 124g n - Xll2g n+ 112 Sn 
Ki = Y V ( lb > 

where X is the average isospin asymmetry (6) of the 
124 Sn- like residue defined as the composition of nucle- 
ons with local densities higher than po/20 and velocities 
larger than 1/2 the beam velocity in the cm. frame. A 
density cut of po/8 is found to give almost same results. 
In ideal case, the value of Ri ranges between 0.05 and 1 
from complete mixing to full transparency. 

A. Effects of momentum-dependent interactions on 
isospin diffusion 

Effects of momentum-dependent interactions on dy- 
namics of heavy-ion collisions can be seen from the time 
evolution of the density distribution [8§|. In Fig. 
we show the density contour p(x, 0, z) in the reaction 
plane at different times for the reaction 124 Sn + 112 Sn at 
E/A = 50 MeV and b = 6 fm calculated with x = — 1 us- 
ing both the MDI and the soft Bertsch-Das Gupta-Kruse 
(SBKD) interactions. It should be noted that the former 



(MDI) interactions are momentum dependent for both 
isoscalar and isovector nuclear potentials while the latter 
(SBKD) interactions do not include any momentum de- 
pendence in either isoscalar or isovector nuclear poten- 
tials though both interactions have similar incompress- 
ibility Kq and the same density dependence in the sym- 
metry energy. The experimental free space JV-JV cross 
sections are used in these calculations. It is seen that 
both interactions give similar time evolution of the col- 
lisions dynamics, namely, projectile-like and target-like 
residues can be clearly separated after about 100 fm/c. 
Detailed examinations indicate, however, that the reac- 
tion system expands more quickly and there are also 
more emitted nucleons in the case of the momentum- 
dependent MDI interaction than that of the momentum- 
independent SBKD interaction. 

In Fig. [5] we show the measured Ri together with the 
predictions from the IBUU04 for the time evolution of 
both Ri and the average central densities calculated with 
x = — 1 using both the MDI and the soft SBKD interac- 
tions. It is seen that isospin diffusion occurs mainly from 
about 30 fm/c to 80 fm/c when average central density 
decreases from about 1.2po to 0.3po- However, the value 
of Ri still changes slightly with time until after about 
120 fm/c when projectile-like and target-like residues are 
well separated as shown in Fig. [7] This is partly due to 
the fact that the isovector potential remains appreciable 
at low density as shown in Fig. [9] where the symmetry 
potential (U n — U p )/25 is shown as a function of nucleon 
momentum (panel (a)) or density (panel (b)) for the MDI 
interaction and as a function of density for the SBKD in- 
teraction (panel (c)). Also, evaluating isospin diffusion 
Ri based on three reaction systems, which have different 
time evolutions for the projectile residue as a result of 
different total energies and numbers of nucleons, further 
contributes to the change of Ri at low density. For the 
two interactions consider here, the main difference be- 
tween the values for Ri appears in the expansion phase 
when densities in the participant region are well below 
po- The experimental data from MSU are seen to be re- 
produced nicely by the MDI interaction with x = — 1, 
while the SBKD interaction with x = — 1 leads to a 
significantly lower value for Ri as the strength of the 
momentum-independent symmetry potential is stronger 
(see Fig. Hp, w hich has been shown to enhance the isospin 
diffusion [92, EM, EM- 

Effects of the symmetry energy on isospin diff usion 
were also studied by varying the parameter x [l09l |. 
Shown in Fig. \W\ is the final saturated value for 1 — R4, 
which measures the degree of isospin diffusion, as a func- 
tion of Kasy for both MDI and SBKD interactions. It is 
obtained by averaging the value of 1 — Ri after 120 fm/c 
with error bars corresponding to its dispersion, whose 
mag nitude is similar to the error band shown in Ref. 
[1081 ] for the theoretical results from the momentum- 
independent BUU model. For the SBKD interaction 
without momentum dependence, the isospin diffusion de- 
creases monotonically (i.e., increasing value for Ri) with 
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FIG. 7: (Color online) Density contour p(x, 0, z) in the reaction plane at different times for the reaction 124 Sn + 112 Sn at 
E I A = 50 MeV and b = 6 fm by using momentum-dependent interaction MDI with x = — 1 (upper panels) and momentum- 
independent interaction SBKD with momentum-independent symmetry potential C/iym 1 ' (p> S, t) (lower panels). Thick solid 
lines represent pn/20 while dashed lines represent po/8. 
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FIG. 8: (Color online) Degree of isospin diffusion as a function 
of time with the MDI and SBKD interactions. Corresponding 
time evolu tions of central density are also shown. Taken from 
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FIG. 9: (Color online) Symmetry potential as a function of 
nucleon momentum (a) or density (b) with the MDI interac- 
tion and SBKD interaction (c). Taken from Ref. |l09l ]. 



increasing strength of K asy as the corresponding isovec- 
tor potential is mostly positive and decreases with in- 
creasing stiffness of E synl (p) in the whole range of consid- 
ered x parameter. The isospin diffusion is reduced when 
the momentum-dependent interaction MDI is used as the 
momentum dependence weakens the strength of symme- 
try potential except for x — —2. As seen in Fig. [9l the 
symmetry potential in the MDI interaction has the small- 
est strength for x — —1 as it is close to zero at k « 1.5 
fm -1 and p/po ~ 0.5, and increases again with further 
hardening of the symmetry energy, e.g., x — —2, when it 
becomes largely negative at all momenta and densities. 
The MDI interaction with x = — 1 thus gives the smallest 



degree of isospin diffusion among the interactions consid- 
ered here and reproduces the MSU data as already shown 
in Fig. El 

The symmetry energy in the MDI interaction with x = 
— 1 can be parameterized as E sym (p) ss 31.6(p/ po) 1 ' 05 . It 
leads to a value of if aS y ~ —550 MeV for the isospin 
dependent part of the isobaric incompressibility of asym- 
metric nuclear matter, which should be compared to the 
published constraint of —566 ±1350 < if a sy < 139 ±1617 
MeV extracted earlier from studying giant monopolc res- 
onances [11| . 
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FIG. 11: (Color online) Time evolution of the isospin diffusion 
Ri using for MDI interactions with different x parameters and 
the i n-me dium nucleon-nucleon cross sections. Taken from 
Ref. Hoi. 



FIG. 10: (Color online) Degree of isospin diffusion as a func- 
tion of -ffasy with the MDI and SBKD interactions. 7 is the 
parameter for fitting the corresponding sy mme try energy with 
E aym (p) = 31.6(p/p ) 7 . Taken from Ref. [Tog| | . 



the experimental data from MSU, which is shown by the 
shaded band. 



B. Effects of in-medium NN cross sections on 
isospin diffusion 

The isospin degree of freedom plays an important role 
in heavy-ion collisions through both the nuclear EOS and 
the nucleon-nucleon (NN) scatterings [3, @]. In partic- 
ular, the transport of isospin asymmetry between two 
colliding nuclei is expected to depend on not only the 
symmetry potential through the density dependence of 
the symmetry energy E sym (p) but also the in-medium 
NN cross sections. The two are related, respectively, to 
the long-range and the short-range parts of the isospin- 
dependent in-medium nuclear effective interactions. For 
instance, the drifting contribution to the isospin trans- 
port in a nearly equilibrium system is proportional to the 
product of the mean relaxation time r np and the isospin 
asymmetric force F np [9l[. While F np is directly related 
to the gradient of the symmetry potential, r np is inversely 
proportional to the neutron-proton (np) scattering cross 
section a np [9l| . Furthermore, the collisional contribu- 
tion to the isospin transport in non-equilibrium system 
is generally expected to be proportional to the np scat- 
tering cross section. In the above study on isospin dif- 
fusion, the free-space NN cross sections were used. In 
this section, we review the effects of in-medium NN cross 
sections on the isospin diffusion in heavy ion collisions. 

Shown in Fig. [11] are the time evolution of Ri re- 
calculated using the in-medium NN cross sections and 
the MDI interaction with four x parameters. It is seen 
that the net isospin transport and the influence of the 
x parameter show up mainly in the expansion phase of 
the reaction after about 40 fm/c and become relatively 
stable after about 80 fm/c. In the late stage of the reac- 
tion, the values of Ri from x = — 1 and come close to 
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FIG. 12: (Color online) Degree of isospin diffusion as a func- 
tion of K aay (po) with the free (filled squares) and in-medium 
(open squares) nucleon-nucleon cross sections. Taken from 
Ref. [i07t . 

For a more meaningful comparison with the experi- 
mental data, the time average of Ri between t = 120 
fm/c and 150 fm/c has been used as in Fig. [8] Shown 
in Fig. [12] is a comparison of the averaged strength of 
isospin transport 1 — Ri obtained with the free and in- 
medium NN cross sections, respectively, as a function of 
the asymmetric part of the isobaric incompressibility of 
nuclear matter at pn, i.e., K asy (po). The error bars were 
drawn to indicate fluctuations and were obta ined from 
the dispersion in the time evolution of Ri jl09l |. First, it 
is interesting to note that with the in-medium NN cross 
sections the strength of isospin transport I—R4 decreases 
monotonically with decreasing value of x. With the free- 
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space NN cross sections, there appears, however, to be a 
minimum at around x = — 1. Moreover, this minimum is 
the point closest to the experimental data. This allowed 
us to extract the value of K asy (po) = —550 ±100 MeV. 
With the in-medium NN cross sections we can now fur- 
ther narrow down the K asy (po) to be about —500 ± 50 
MeV. The latter is consistent with that extracted from 
studying the isospin dependence of giant resonan ces o f 
112 Sn to 124 Sn isotopes by Fujiwara et al at Osaka |l59| . 
Also shown in the figure are the 7 values used in fitting 
the symmetry energy with E sym {p) = 31.6(p/po) 7 - The 
results with the in-medium NN cross sections constrain 
the 7 parameter to be between 0.69 and 1.05. The lower 
value is clo se to wha t is extracted from studying giant 
resonances jl60l . Il6l| . The value of 7 = 1.05 extracted 
earlier using the free-space NN cross sections sets an up- 
per limit. 

The difference in 1 — Ri obtained with the free-space 
and the in-medium NN cross sections is small for x = 1 
and x = 0, but becomes large for x = — 1 and x = — 2. 
The increasing effect of the in-medium NN cross sections 
with decreasing K asy (po) or x parameter can be under- 
stood from consideration of the contributions from the 
symmetry potential and the np scatterings. As we have 
mentioned above, both contributions to the isospin trans- 
port depend on the np scattering cross section a np . While 
the collisional contribution is proportional to the np scat- 
tering cross section <7 npi the mean-field contribution is 
proportional to the product of the isospin asymmetric 
force F np and the inverse of a np . The overall effect of the 
in-medium NN cross sections on isospin transport is a re- 
sult of a complicated combination of both the mean field 
and the NN scatterings. Generally speaking, the symme- 
try potential effects on the isospin transport will become 
weaker when the NN cross sections are larger while the 
symmetry potential effects will show up more clearly if 
smaller NN cross sections are used. This feature can be 
seen from Fig.[9l For x = 1 and x = 0, the symmetry po- 
tential and its gradient with respect to density, as shown 
in Fig. [21 are large at low densities where the majority of 
net isospin transport occurs. The F np factor makes the 
contribution due to the mean field to dominate over that 
due to the collisions. Therefore, the reduced in-medium 
<j np leads to about the same or a slightly higher isospin 
transport. As the x parameter decreases to x = — 1 and 
x = —2, however, the symmetry potential decreases and 
its density gradient can be even negative at low densi- 
ties. In these cases, either the collisional contribution 
dominates or the mean-field contribution becomes neg- 
ative. Therefore, the reduced in-medium np scattering 
cross section a np leads to a weaker isospin transport com- 
pared with the case with the free-space NN cross sections. 

Based on an isospin- and momentum-dependent 
IBUU04 transport model with free-space experimental 
NN cross sections, comparing the theoretical results with 
the experimental data has allowed us to extract a nuclear 
symmetry energy of E sym (p) ss 31.6(p/po) ■ Including 
also medium-dependent NN cross sections, which are im- 
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FIG. 13: (Color online) Degree of the isospin diffusion 1 — 
Ri as a function of L using the MDI interaction with x = 
—2, —1, 0, a nd 1 . The shaded band indicates the data from 
NSCL/MSU |l08t |. The solid square with error bar represents 
L = 88 ± 25 MeV. Taken from Ref. [Til. 



portant for isospin-dependent observables [l07l Il62| , the 
isospin diffusion data leads to an even softer nuclear sym- 
metry energy of E sym (p) w 31.6(p/po) 7 with 7 w 0.69 
[l07j . 

Since the slope parameter L of the nuclear symmetry 
energy gives an important constraint on the density de- 
pendence of the nuclear symmetry energy, in particular, 
it has been shown that the slope parameter L is related 
to the neutron skin thickness of heavy nuclei, it is in- 
teresting to see how the isospin diffusion data constrain 
the value of L. In Fig. [IBJ we show the results from the 
IBUU04 transport model with in-medium NN cross sec- 
tions, that are consistent with the mean-field potential 
obtained with the MDI interactions used in the model, 
for the degree of the isospin diffusion 1 — Ri as a function 
of L. The sh aded band in Fig. [TBI indicates the data from 
NSCL/MSU [l08t . It is seen that the strength of isospin 
diffusion 1 — Ri decreases monotonically with decreasing 
value of x or increasing value of L. This is expected as 
the parameter L reflects the difference in the pressures 
on neutrons and protons. From comparison of the theo- 
retical results with the data, we can clearly exclude the 
MDI interaction with x = 1 and x = — 2 as they give ei- 
ther too large or too small a value for 1 — Ri compared to 
that of data. The range of x or L values that give values 
of 1 — Ri falling within the band of experimental values 
could in principle be determined in our model by detailed 
calculations. Instead, we determine this schematically by 
using the results from the four x values. For the centroid 
value of L, it is obtained from the interception of the line 
connecting the theoretical results at x = — 1 and with 
the central value of 1 — R4 data in Fig. [131 i.e., L = 88 
MeV. The upper limit of L = 113 MeV is estimated from 
the interception of the line connecting the upper error 
bars of the theoretical results at x = — 1 and —2 with 
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the lower limit of the data band of 1 — R4. Similarly, the 
lower limit of L = 65 MeV is estimated from the inter- 
ception of the line connecting the lower error bars of the 
theoretical results at x = and —1 with the upper limit 
of the data band of 1 — R4. This leads to an extracted 
value of L — 88 ± 25 MeV as shown by the solid square 
with error bar in Fig. 1131 



V. CONSTRAINING THE SKYRME 
EFFECTIVE INTERACTIONS AND THE 
NEUTRON SKIN THICKNESS OF NUCLEI 
USING ISOSPIN DIFFUSION DATA FROM 
HEAVY ION COLLISIONS 

Information on the density dependence of the nuclear 
symmetry energy can in principle also be obtained from 
the thickness of the neutron skin in heavy nuclei as the 
latter is strongly correlated with the slope parameter L 
of t he nuclear symmet r y en ergy at saturation density 
O, EM, EH ml EH [H^. Because of the large un- 
certainties in measured neutron skin thickness of heavy 
nuclei, this has not been possible. Instead, studies have 
been carried out to use the extracted nuclear symmetry 
energy from the isospin diffusion data t o constra in the 
neutron skin thickness of heavy nuclei [l07l . IllOj ], Us- 
ing the Hartree-Fock approximation with parameters fit- 
ted to the phenomenological EOS that was used in the 
IBUU04 transport model to describe the isospin diffusion 
data from the NSCL/MSU, it was f ound that a neutron 
skin thickness of less than 0.15 fm [IH |TT3| for 208 Pb 
was incompatible with the isospin diffusion data. 

In the following, we study more systematically the cor- 
relation between the density dependence of the nuclear 
symmetry energy and the thickness of the neutron skin in 
a number of nuclei within the framework of the Skyrme 
Hartree-Fock model. Using the extracted values of L 
from the isospin diffusion data in heavy- ion collisions, we 
obtain stringent constraints on the neutron skin thick- 
ness of the nuclei 208 Pb, 132 Sn, and 124 Sn. The extracted 
value of L also limits the allowed parameter sets for the 
Skyrme interaction. 



A. Constraining the Skyrme effective interactions 
using isospin diffusion data 

In the standard Skyrme Hartree-Fock model, the in- 
teraction is taken to have a z ero-range , density- and 
momentum-dependent form (oil . I168L Il69l ll70L Il7l| . i.e., 

Ui 2 (R,r) = i (l + x P a )6(r) 

+ ± : t 3 (l + x 3 P„)p' T (R.)6{r) 
6 

+ ^t 1 {\ + x 1 P (y ){K' 2 8{Y) + 5{v)K 2 ) 

+ t 2 {l + x 2 P a )K -S(t)K 
+ iW K' -aCrJKoi+o*) x K], 



- Skyrme-Hartree-Fock 
with 21 parameter sets 
MDI interaction with x= 
MDI interaction with x= -1 




FIG. 14: (Color online) Density dependence of the nuclear 
symmetry energy E ayln (p) for 21 sets of Skyrme interaction 
parameters. The results from the MDI interaction with x = 
— 1 (open s quar es) and (solid squares) are also shown. Taken 
from Ref. [ml ]. 



with r = i - ! — r 2 and R = (r x + r 2 )/2. In the above, 
the relative momentum operators K = (Vi — V 2 )/2i and 
K' = — (Vi — V 2 )/2i act on the wave function on the right 
and left, respectively. The quantities P a and <Ji denote, 
respectively, the spin exchange operator and Pauli spin 
matrices. The a, to — t 3 , xq — x 3 , and Wo are Skyrme 
interaction parameters that are chosen to fit the bind- 
ing energies and charge radii of a large number of nuclei 
in the periodic table. For infinite nuclear matter, the 
symmetry en ergy from the Skyrme interaction can be 
expressed as [l70L Il7l| 



lh 2 Att 2 ^ 73 



Es ^ ip) = 3 2^ {— 

- h (2x + l)p 

2 \ 2 /3 



2/3 



48 



t 3 (2x 3 + l)p° 



1 /3tH , r 
+ (4 + 5x 2 )i 2 ]p 5 / 3 . 



(17) 



(16) 



The coefficient of the S 4 term in Eq. ([I]) can also be ob- 
tained analytically and has been shown to be very small 
over a large range of nuclear density (< 0.80 fm -3 ) and 
isospin asymmetry. The parabolic law of Eq. |T]) wi thou t 
the r5 4 and higher-order terms in d is thus justified [l70l |. 

Fig. [T3] displays the density dependence of E syin (p) 
for 21 sets of Skyrme interaction parameters, i.e., SKM, 
SKM*, RATP, SI, SI I, SI II, SIV, SV, SVI, E, E a , 
G cr , i?cr, Z, Z (j , Z*, T, T3, SkX, SkXce, and SkXm. 
The values of the parameters i n these S kyrme interac- 
tions can be found in Refs. [6li I168L Il69| |. For compari- 
son, we also show in Fig. [T3]results from the phenomeno- 
logical MDI interactions with x = — 1 (open squares) and 
(solid squares). As we have discussed above, from com- 
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paring the isospin diffusion data from NSCL/MSU using 
the IBUU04 with in-medium NN cross sections, these in- 
teractions are recently shown to give, respectively, the 
upper and lowe r bounds for the stiffness of the symme- 
try energy 107(. It is seen from Fig. [14] that the den- 
sity dependence of the symmetry energy varies drasti- 
cally among different interactions. Although the values 
of E sym (po) are all in the range of 26-35 MeV, the val- 
ues of L and K sym are in the range of —50-100 MeV and 
-700-50 MeV, respectively. 

The extracted value of L = 88 ± 25 MeV gives a rather 
stringent constraint on the density dependence of the nu- 
clear symmetry energy and thus puts strong constraints 
on the nuclear effective interactions as well. For the 
Skyrme effective interactions shown in Fig. I14[ for in- 
stance, all of those lie beyond x = and x = — 1 in 
the sub-saturation region are not consistent with the ex- 
tracted value of L. Actually, we note that only 4 sets of 
Skyrme interactions, i.e., SIV, SV, G CT , and R CT , in the 21 
sets of Skyrme interactions considered here have nuclear 
symmetry energies that are consistent with the extracted 
L value. 



B. Constraining the neutron skin thickness of 
nuclei using isospin diffusion data 



The neutron skin thickness S of a nucleus is defined as 
the difference between the root-mean-square radii \J (r n ) 
of neutrons and U (r p ) of protons, i.e., 



(18) 



It has been known that S is sensitive to the density de- 
pendence of the nuclear symmetry energy, particularly 
the slope paramet e r L at the norm al nuclear matter den- 
sity [lfl, UM, EH, EH, EH EH]- Using above 21 sets 
of Skyrme interaction parameters, we have evaluated the 
neutron skin thickness of several nuclei. In Figs. fTBT a). 
(b) and (c), we show, respectively, the correlations be- 
tween the neutron skin thickness of 208 Pb with L, K sym , 
and E sym (po). It is seen from Fig. [ToTa) that there exists 
an approximate linear correlation between S and L. The 
correlations of S with K sym and E sym (po) are less strong 
and even exhibit some irregular behavior. The solid line 
in Fig. [TSTa') is a linear fit to the correlation between S 
and L and is given by the following expression: 

5( 208 Pb) = (0.1066 ±0.0019) 

+ (0.00133 ± 3.76 x 10~ 5 ) x L, (19) 



L = (-78.5 ±3.2) 

+ (740.4 ±20.9) x 5( 208 Pb), 



(20) 



where the units of L and S are MeV and fm, respectively. 
Therefore, if the value for either 5'( 208 Pb) or L is known, 
the value for the other can be determined. 
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FIG. 15: (Color online) Neutron skin thickness S of 208 Pb 
as a function of (a) L, (b) K sym , and (c) E syln (po) for 21 
sets of Skyrme interaction parameters. The line in panel (a) 
represents a linear fit. Taken from Ref. [l!2j |. 



It is of interest to see if there are also correlations be- 
tween the neutron skin thickness of other neutron-rich 
nuclei and the nuclear symmetry energy. Fig. [TH] shows 
the same correlations as in Fig. [15] but for the neutron- 
rich nuclei 132 Sn, 124 Sn, and 48 Ca. For the heavy 132 Sn 
and 124 Sn, we obtain a similar conclusion as for 208 Pb, 
namely, S exhibits an approximate linear correlation with 
L but weaker correlations with -Ksym and -E S ym(po)- For 
the lighter 48 Ca, on the other hand, all the correlations 
become weaker than those of heavier nuclei. Therefore, 
the neutron skin thickness of heavy nuclei is better cor- 
related with the density dependence of the nuclear sym- 
metry energy. As in Eqs. (fT!?]) and (|2"0|) , a linear fit to 
the correlation between S and L can also be obtained for 
132 Sn and 124 Sn, and the corresponding expressions are 



x L. 



and 



S'( 132 Sn) = (0.1694 ± 0.0025) 

+ (0.0014 ±5.12 x 10~ 5 ) 

L = (-117.1 + 5.4) 

+ (695.1 + 25.3) x S-( 132 Sn), 



S'( 124 Sn) = (0.1255 + 0.0020) 

+ (0.0011 + 4.05 x 10~ 5 ) x L. 

L = (-110.1 + 5.2) 

+ (882.6 + 32.3) x S , ( 124 Sn), 



(21) 



(22) 



(23) 



(24) 



Similar linear relations between S and L are also ex- 
pected for other heavy nuclei. This i s not surprising as 
detailed discussions in Refs. 0, EH, EH, ITdIIIdI fl67j 
have shown that the thickness of the neutron skin in 
heavy nuclei is determined by the pressure difference be- 
tween neutrons and protons, which is proportional to the 
parameter L. 

To give a quantitative estimate of above discussed cor- 
relations, we define the following linear correlation coef- 
ficient C;: 



Ci = y/1 - q/t, 



(25) 
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FIG. 16: (Color online) Same as Fig. 2 but for nuclei 132 Sn 
(Solid squares), 1 24 Sn (Open squares) and 48 Ca (Triangles). 
Taken from Ref. [Hi]. 



TABLE II: Linear correlation coefficients C; of 5* with L, K syln 
and E Byrn ( Po ) for 208 Pb, 132 Sn, 124 Sn, and 48 Ca from 21 sets 
of Skyrme interaction parameters. Taken from Ref. [I12j |. 



a (%) 


208 p b 


132 Sn 


124 Sn 


48 Ca 


S-L 


99.25 


98.76 


98.75 


93.66 


S-Ksym 


92.26 


92.06 


92.22 


86.99 


S-Esym 


87.89 


85.74 


85.77 


81.01 



where 

n 

i = ^[y,-(A + i?^)] 2 , (26) 

n n 

i=i i=i 

In the above, A and B are the linear regression coeffi- 
cients, (xi, yi) are the sample points, and n is the num- 
ber of sample points. The linear correlation coefficient 
Ci measures the degree of linear correlation, and Ci — 1 
corresponds to an ideal linear correlation. Table [Til gives 
the linear correlation coefficient C; for the correlation of 
S with L, K sym and E sym ( Po ) for 208 Pb, 132 Sn, 124 Sn, 
and 48 Ca shown in Figs. [TBI and [TBI for different Skyrme 
interactions. It is seen that these correlations become 
weaker with decreasing nucleus mass, and a strong lin- 
ear correlation only exists between the S and L for the 
heavier nuclei 208 Pb, 132 Sn, and 124 Sn. Therefore, the 
neutron skin thickness of these nuclei can be extracted 
once the slope parameter L of the nuclear symmetry en- 
ergy at saturation density is known. 

The extracted L value from isospin diffusion data al- 
lows us to determine from Eqs. (fT5)) . ([2"T)l . and ([25)1 . 
respectively, a neutron skin thickness of 0.22 ± 0.04 fm 
for 208 Pb, 0.29 ± 0.04 fm for 132 Sn, and 0.22 ± 0.04 fm 
for 124 Sn. Experimentally, great efforts were devoted to 
mea sure the t hickness of the neutron skin in heavy nu- 
clei [l72|, Il73j , and a recent review can be found in Ref. 

The data for the neutron skin thickness of 208 Pb 
indicate a large uncertainty, i.e., 0.1-0.28 fm. Our re- 
sults for the neutron skin thickness of 208 Pb are thus 



consistent with present data but give a much stronger 
constraint. A large uncertainty is also found experimen- 
tally in the neutron skin thickness of 124 Sn, i.e., its value 
varies from 0.1 fm to 0.3 fm depending on the exper- 
imental method. The proposed experiment of parity- 
violating electron scattering from 208 Pb at the Jefferson 
Laboratory is expected to give another independent and 
more accurate measurement of its neutron skin thickness 
(within 0.05 fm), thus providing improved constraints on 
the densi ty dependence of the nuclear symmetry energy 

[MUzl- 

Recently, an accurately calibrated relativistic 
parametrization based on the relativistic mean-field 
theory has been introdu ced to study the neutron skin 
thickness of finite nuclei jll3| . This parametrization can 
describe simultaneously the ground state properties of 
finite nuclei and their monopole and dipole resonances. 
Using this parametrization, the authors predicted a 
neutron skin thickness of .21 fm in 208 Pb, 0.27 fm in 
132 Sn, and 0.19 fm in 124 Sn [HEEIl- These predictions 
are in surprisingly good agreement with our results 
constrained by the isospin diffusion data in heavy-ion 
collisions. 

In addition, the neutron skin thickness of the nucleus 
90 Zr has recently been determined to be 0.07 ± 0.04 fm 
from the model-independent spin-dipole sum rule value 
measured from the charge-exchange spin-dipole excita- 
tions [177} . This value is reproduced by the symmetry 
energy with L = 88 ± 25 MeV extracted from the isospin 
diffusion data in heavy-ion collisions, which predicts a 
neutron skin thickness of 0.088 ± 0.04 fm for 90 Zr. 



VI. PROBING THE HIGH DENSITY 
BEHAVIOR OF THE NUCLEAR SYMMETRY 

ENERGY IN HEAVY-ION REACTIONS 
INDUCED BY HIGH ENERGY RADIOACTIVE 
BEAMS 

Although significant progress has been made in the 
determination of the density dependence of the nuclear 
symmetry energy at sub-normal densities, the high den- 
sity behavior of the the nuclear symmetry energy is still 
poorly known. Fortunately, heavy-ion reactions, espe- 
cially those induced by high energy radioactive beams to 
be available at high energy radioactive beam facilities, 
provide a unique opportunity to pin down the high den- 
sity behavior of the symmetry energy. In this section, 
we illustrate via transport model simulations several ex- 
perimental observables which are sensitive to the high 
density behavior of the symmetry energy. 



A. Isospin asymmetry of dense matter formed in 
high energy heavy-ion reactions 

To see the maximum baryon density and isospin asym- 
metry that can be achieved in central heavy ion colli- 
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FIG. 17: (Color online) Central baryon density (upper win- 
dow) and isospin asymmetry (lower window) of high density 
region for the reaction of 132 Sn -|- 124 Sn at a beam energy of 
400 MeV/nucleon and an impact parameter of 1 fm. Taken 
from Ref. 0. 



sions induced by high energy radioactive beams in fu- 
ture radioactive beam facilities, we show in Fig. [17] the 
central baryon density (upper window) and the average 
(n/p)p> Po ratio (lower window) of all regions with baryon 
densities higher than pg in the reaction of 132 Sn+ 124 Sn 
at a beam energy of 400 MeV/nucleon and an impact 
parameter of 1 fm. It is seen that the maximum baryon 
density is about 2 times normal nuclear matter density. 
Moreover, the compression is rather insensitive to the 
symmetry energy because the latter is relatively small 
compared to the EOS of symmetric matter around this 
density. The high density phase lasts for about 15 fm/c 
from 5 to 20 fm/c for this reaction. The isospin asym- 
metry of the high density region is, however, sensitive to 
the symmetry energy. The soft (e.g., x = 1) symmetry 
energy leads to a significantly higher value of (n/p) p > Po 
than the stiff one (e.g., x = —2). This is consistent 
with the well-known isospin fractionation phenomenon 
that it is energetically more favorable to have a higher 
isospin asymmetry 5 in the high density region for a 
softer symmetry energy functional E sym (p) as a result 
of the E sym (p)5 2 term in the EOS of asymmetric nuclear 
matter. Since the symmetry energy changes from be- 
ing soft to stiff when the parameter x varies from 1 to 
—2, the value of (n/p) p > Po becomes lower in the supra- 
normal density region as the parameter x changes from 
1 to —2. Because of neutron-skins of the colliding nu- 
clei, especially that of the projectile 132 Sn, the initial 
value of the quantity (n/p) p > Po , which is about 1.4, is 
less than the average n/p ratio of 1.56 of the reaction 
system. Also, in neutron-rich nuclei, the n/p ratio on 
the low-density surface is much higher than that in their 



interior. The dense matter region in heavy ion collisions 
can thus become either neutron-richer or neutron-poorer 
with respect to the initial state depending on the sym- 
metry energy functional E syin (p) used. 



B. Isospin fractionation and n-p differential flow 

The degree of isospin equilibration or translucency 
in heavy ion collisions can be measured by the rapid- 
ity distribution of nucleon isospin asymmetry <5f rcc = 
(N n - N p )/(N n + N p ), where N n and N p are multi- 
plicities of free neutrons and protons, respectively [90| . 
Although it might be difficult to measure directly 5{ Tee 
because it requires the detection of neutrons, similar in- 
formation can be extracted from ratios of light clusters, 
such as, t/ 3 He, as demonstrated recently within a coa- 
lescence model [86, 89]. Shown in Fig. [TS] are rapidity 
distributions of Sf Toe with (upper window) and without 
(lower window) the Coulomb potential. It is interest- 
ing to see that the 5f lee at midrapidity is particularly 
sensitive to the symmetry energy. As the parameter x 
increases from —2 to 1 the i5f ree at midrapidity decreases 
by about a factor of 2. Moreover, the forward-backward 
asymmetric rapidity distributions of <5f roc with all four x 
parameters indicates the apparent nuclear translucency 
during the reaction [95| . 
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FIG. 18: (Color online) Isospin asymmetry of free nucleons 
with and without the Coulomb force for different symmetry 
energies. Taken from Ref. [95j . 

Another observable that is sensitive to the high density 
behavior of the symmetry energy is the neutron-proton 
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differential flow first introduced in Ref. 8z 



n—p 



(28) 



i=l 



where Wi = 1(— 1) for neutrons (protons) and N(y) is 
the total number of free nucleons at rapidity y. Since 
the differential flow depends on the symmetry potential 
through the latter's effects on the isospin fractionation 
and the collective flow, it has the advantage of maximiz- 
ing the effects of the symmetry potential while minimiz- 
ing those of the isoscalar potential. Shown in Fig. [19] is 
the n-p differential flow for the reaction of 132 Sn+ 124 Sn 
at a beam energy of 400 MeV/nucleon and an impact 
parameter of 5 fm. Effects of the symmetry energy are 
clearly revealed by changing the parameter x. 
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FIG. 20: (Color online) The vr" and ^ 
the x parameter. Taken from Ref. 



yields as functions of 



rich for the softer symmetry energy as a result of isospin 
fractionation [9~jj . the ir~ multiplicity is thus more sensi- 
tive to the isospin asymmetry of the reaction system and 
the symmetry energy. However, it is well known that the 
pion yield is also sensitive to the symmetric part of the 
nuclear EOS, and it is thus hard to extract reliable in- 
formation about the symmetry energy from the ir~ yield 
alone. The ir~/ir + ratio is, on the other hand, a bet- 
ter probe as this ratio is sensitive only to the diffe renc e 
in the chemical potentials for neutrons and protons [l79j | . 
As well demonstrated in Fig. [3T1 the 1 7r+ ratio is quite 
sensitive to the symmetry energy, especially at low trans- 
verse momenta, and can be used to probe the high density 
behavior of nuclear symmetry energy E sym (p). 



FIG. 19: (Color online) Neutron-proton differ entia l flow for 
different symmetry energies. Taken from Ref. [l03| ]. 



C. Pion yields and tt /n + ratio 

At beam energy of about 400 MeV/nucleon, pion pro- 
duction becomes non-negligible and may also carry inter- 
esting infor mati on about the EOS of dense neutron-rich 
matter [iJ, Il78| . In Fig. [201 we show the ir~ and tt + 
yields as functions of the x parameter in the MDI inter- 
action. It is seen that when the x parameter is changed 
from -2 to 1, the 7r~ multiplicity increases by about 20%, 
although the ir + multiplicity remains about the same. 
The 7T _ multiplicity is thus more sensitively to the sym- 
metry energy than that of ir + . Also, the multiplicity of 
7r~ is about 2 to 3 times that of 7r + . This is because tt~ 
mesons are mostly produced from neutron-neutron colli- 
sions, which are more frequent in collisions of neutron- 
rich nuclei. Since the high density region is more neutron 



3.6 i 1 1 1 1 r 




1.6 ' 1 ' ' 1 1 1 1 >— 

50 100 150 200 250 



p t cm (MeV/c) 

FIG. 21: (Color online) The ix~ ratio as a function of 
transverse momentum. Taken from Ref. [§3|. 
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D. Double n/p and it ratio 

Because of the facts that the symmetry potentials have 
opposite signs for neutrons and protons and are also 
generally smaller compared to the isoscalar potential at 
the same density, most observables proposed so far for 
studying the density dependence of the nuclear symme- 
try energy employ differences or ratios of isospin multi- 
plets of baryons, mirror nuclei and mesons, such as, the 
neutron/proton ratio of emitted nucleons [681 ] . neutron- 
proton differential flow 1821 . neutron- pro ton corre l ation 
function Jp, t/ 3 Re \g&W&, ^"At ^ ^ 
Xr/£+ [9J and K°/K+ ratios [l02j], etc. Among these 
observables, the ratio of emitted neutrons to protons 
has probably the highest sensitivity to the symmetry en- 
ergy as the symmetry potential acts directly on nucleons 
and emitted nucleons are also rather abundant in typi- 
cal heavy-ion reactions. However, it is very challenging 
to measure some of these observables, especially those 
involving neutrons. The measurement of neutrons, par- 
ticularly the low energy ones, always suffers from low 
detection efficiencies even for the most advanced neu- 
tron detectors. Therefore, observables involving neutrons 
normally have large systematic errors. Moreover, for 
essentially all of these observables, the Coulomb force 
on charged particles plays an important role and some- 
times competes strongly with the symmetry potential. 
One has to disentangle carefully effects of the symme- 
try potential from those due to the Coulomb potential. 
It is thus very desirable to find experimental observables 
which can reduce the influence of both the Coulomb force 
and the systematic errors associated with neutrons. A 
possible candidate for such an observable is the dou- 
ble ratios of emitted neutrons and protons taken from 
two reaction systems using four isotopes of the same el- 
ement, namely, the neutron/proton ratio in the neutron- 
rich system over that in the more sym metric system, as 
recently proposed by Lynch et al. [l80j ] . They have actu- 
ally demonstrated the feasibility of measuring the double 
neutron/proton ratios in central reactions of 124 Sn+ 124 Sn 
and 112 Sn+ 112 Sn at a beam energy of 50 MeV/nucleon 
at th e National Superconducting Cyclotron Laboratory 

pot 

Both the double neutron/proton ratio and the dou- 
ble 7T-/7T+ ratio in 132 Sn+ 124 Sn and 112 Sn+ 112 Sn reac- 
tions at 400 MeV/nucleon have been studied in the IBUU 
model in order to demonstrate the effect of symmetry en- 
ergy at high density. It was found that these ratios have 
about the same sensitivity to the density dependence of 
symmetry energy as the corresponding single ratio in the 
respective neutron-rich system involved. Given the ad- 
vantages of measuring the double ratios over the single 
ones, the study of double ratios will be more useful for 
further constraining the symmetry energy of neutron-rich 
matter. Furthermore, the systematic errors associated 
with transport model calculations are mostly related to 
the uncertainties in the in-medium NN cross sections, 
techniques of treating collisions, sizes of the lattices in 



calculating the phase space distributions, techniques in 
handling the Pauli blocking, etc. Since the double ratio is 
a relative observable from two similar reaction systems, 
these systematic errors are expected to be reduced. 
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FIG. 22: (Color online) The double neutron/proton ratio 
of free nucleons in the reaction of 132 Sn + 124 Sn and 112 
Sn+ 112 Sn at 400 MeV/nucleon and an impact parameter of 1 
fm (left wi ndow ) and 5 fm (right window), respectively. Taken 
from Ref. p^ . 

In Fig. [22l we show the double neutron/proton ratios 
from the reactions of 132 Sn+ 124 Sn and 112 Sn+ 112 Sn at 
a beam energy of 400 MeV/nucleon and an impact pa- 
rameter of 1 fm (left window) and 5 fm (right window) 
predicted by the IBUU model using the MDI interaction 
with x = and x = — 1, which are consistent with the 
symmetry energy used for sub-saturation densities. At 
both impact parameters, effects of the symmetry energy 
are about 5% — 10% changing from the case with x = to 
x = — 1. One notices here that the low energy nucleons 
have the largest sensitivity to the variation of the sym- 
metry energy for such high energy heavy-ion collisions. 
In fact, the neutron/proton ratio of midrapidity nucleons 
which have gone through the high density phase of the 
reaction are known to be most sensitive to the symmetry 
energy [95| . Compared to the results at the beam energy 
of 50 MeV/nucleon [95j], it is interesting to see a clear 
turnover in the dependence of the double neutron/proton 
ratio on the x parameter, namely the double ratio is lower 
at 50 MeV/nucleon but higher at 400 MeV/nucleon with 
x = — 1 than that with x = 0. The maximum density 
reached at the beam energy of 50 and 400 MeV/nucleon 
is about 1.2po and 2po [951 ], respectively. The turnover 
clearly indicates that the double neutron/proton ratio 
reflects closely the density dependence of the symmetry 
energy as shown in Fig.O This observation also indicates 
that systematic studies of the double neutron/proton ra- 
tio over a broad beam energy range will be important for 
mapping out the density dependence of the symmetry 
energy. 

As shown before from both the total yields and the mo- 
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FIG. 23: (Color online) Kinetic energy distribution of the 
single 7T-/7T+ ratio for 132 Sn+ 124 Sn and 112 Sn+ 112 Sn at 
a beam energy of 400 MeV/nucleon and an impact param- 
eter of b — 1 fin with the stiff (x — —1) and soft (x = 0) 
symmetry energies. The dash-dotted line is the single ir~ 

ratio obtained by turning off the Coulom b po tentials in the 
132 Sn+ i24 gn reaction> Xaken from Re f. fig^l. 



about 2.4. This value agrees with the predicted value 
of (57V 2 + NZ)/(5Z 2 + NZ) w (N/Z) 2 » 2.43, where 
N and Z are the total neutron and proton numbe rs in 
the participant region, from the A resonance model 181 1 
for central 132 Sn+ 124 Sn collisions. This is not surpris- 
ing as at 400 MeV/nucleon, pio ns a re almost exclusively 
produced via the A resonances |l82j | . Comparison of cal- 
culated results with and without the Coulomb potentials 
thus indicates clearly that the peak observed in the sin- 
gle ir~ /ir + ratio is due to the Coulomb effects. Although 
the Coulomb potential distorts the spectra of pions, the 
effect of symmetry potential is still seen in the result- 
ing 7r~/7r + , particularly near its peak value where pions 
have relatively low kinetic energies. These pions are pro- 
duced in the high density nucleonic matter (about 2po) 
through the A resonances and experience many rescat- 
terings with nucleons at both high and low densities as 
well as the Coulomb potential from protons at different 
densities. Since the soft (x = 0) and stiff (x — — 1) sym- 
metry energies have slight difference at low densities and 
the large difference appears at high densities (about 2p ) 
as shown in Fig. [21 the observed symmetry energy effects 
on the energy dependence of the tt~ /ir + ratio thus mainly 
reflect (though not completely) informati on on the high 
density behavior of the symmetry energy [ 1 lj ] . 



mentum spectra, the tt~ /ir + ratio is a promising probe 
of th e symmetry energy at high densities [H, [H, [93, llOOl 
Il0l| . In Fig. [231 we show again the kinetic energy dis- 
tribution of the single 7r _ /7r + ratio for the reactions of 
i32 Sn+ i24 Sn and ii2 Sn+ ii2 Sn at a beam energy f 400 

MeV/nucleon and an impact parameter of b = 1 fm with 
the stiff (x = —1) or soft (x = 0) symmetry energy, ob- 
tained with 12000 IBUU events for each reaction. It is 
seen that the overall magnitude of it~ /tt + ratio is larger 
for the neutron-rich system 132 Sn+ 124 Sn than for the 
neutron-deficient system 112 Sn+ 112 Sn as expected. Also, 
the soft symmetry energy (x = 0) leads to a larger single 
7r~/7r + ratio than the stiff one (x = —1). This is mainly 
because the high density region where most pions are 
produced are more neutron-rich with the softer symme- 
try energy as a result of the isospin fractionation 0, [94[ . 
Furthermore, the single 7r~ /n + ratio is more sensitive to 
the symmetry energy in the reaction 132 Sn+ 124 Sn than in 
the reaction 112 Sn+ 112 Sn as a result of the larger isospin 
asymmetry in the more neutron-rich system. 

Fig. [23] shows that the single ir~ /ir + ratio exhibits 
a peak at a pion kinetic energy of about 45 MeV in 
all cases considered here. The origin of this peak can 
be understood from the single tt~/t: + ratios in both 
reactions from turning off the Coulomb potentials for 
all charged particles. As an example, shown in Fig. 
[23l with the dash-dotted line is the single ir~ /ir + ra- 
tio obtained by turning off the Coulomb potentials in 
the 132 Sn+ 124 Sn reaction. It is seen that the single 
7r~ /tt + ratio now becomes approximately a constant of 
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FIG. 24: (Color online) Kinetic energy dependence of the 
double 7T-/7T+ ratio of 132 Sn+ 124 Sn over 112 Sn+ 112 Sn at a 
beam energy of 400 MeV/nucleon and an impact parameter 
b = 1 f m with the stiff (x = — 1) and soft (x — 0) symmetry 
energies. Taken from Ref. [I84j |. 

In order to reduce the systematic errors related to the 
symmetry energy effect on the tt~ /ir + ratio, it is more 
useful to study the double it~ /tt + ratio in the reactions 
of 132 Sn+ 124 Sn and 112 Sn+ 112 Sn as for the double pro- 
ton/neutron ratio. Fig. EMI shows the double ir~ /ir + ratio 
for these two reactions. It is seen that the kinetic energy 
dependence of the double tt" /ir + ratio is rather differ- 
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ent for the stiff (x = —1) and soft (x — 0) symmetry 
energies. While the double ir~ /ir + ratio is quite flat for 
x = 0, it displays a concave structure for x = — 1 around 
the Coulomb peak. These different behaviors can be un- 
derstood from the corresponding single ir~ /ir + ratios in 
the two reactions shown in Fig. [23l It is reassuring to see 
that around the Coulomb peak the double it~ /tt + ratio 
is still sensitive to the symmetry energy. Compared with 
the single 7r~ /tt + ratio, the kinetic energy dependence of 
the double tt~ /tt + ratio becomes, however, weaker. This 
is because effects of the Coulomb potentials are reduced 
in the double tt^/tt + ratio. We note that the double 
7r~ /ir + ratio displays an opposite symmetry energy de- 
pendence compared with the double n/p ratio for free 
nucleons shown in Fig. [22l This is understandable since 
the soft symmetry energy leads to a more neutron-rich 
dense matter in heavy-ion collisions induced by neutron- 
rich nuclei and thus a smaller n/p ratio for free nucleons 
due to the charge conservation. On the other hand, more 
7r~'s would be produced due to more neutron-neutron in- 
elastic scatterings in the more neutron-rich matter. 



E. Double neutron-proton differential transverse 
flow 



The neutron-proton differential transverse flow defined 
in Eq. (|28|) can be further expressed as [H, HH, l94j | 



reduced to 



K- P (y) 



i 



N(y) 



p*^ 1 



<»> 



where N(y), N n (y) and N p (y) are the numbers of free 
nucleons, neutrons, and protons, respectively, at rapidity 
y\ Piiy) is the transverse momentum of a free nucleon 
at rapidity y; Wi = 1 (—1) for neutrons (protons); and 
(Pn(y)) an d (Pp{y)) are > respectively, the average trans- 
verse momenta of neutrons and protons at rapidity y. Eq. 
([29]) shows that the neutron-proton differential transverse 
flow depends not only on the proton and neutron trans- 
verse momenta but also on their relative multiplicities, 
i.e., the isospin fractionation. This can be more clearly 
seen from the following two special cases. If neutrons 
and protons have the same average transverse momen- 
tum in the reaction plane but different multiplicities in 
each rapidity bin, i.e., (p x n {y)) = (p p (y)} = (p x (y)}, and 
N n {y) ^ N p (y), then Eq. {25} is reduced to 



K- P (y) 



N n (y) - N p (y) 
N(y) 



{p x (y))=S(y)-( P x (y)), (30) 



reflecting the effects of isospin fractionation. On the 
other hand, if neutrons and protons have the same mul- 
tiplicity but different average transverse momenta, i.e., 
N n (y) = N p (y) but (p^y)) £ {p x (y)), then Eq. (25) is 



K- P (y) = ~((p x n (y)) - (p x (y))) 



(31) 



and reflects directly the difference between the neutron 
and proton transverse flows. Since a stiff er symmetry 
potential is expected to lead to a higher isospin fraction- 
ation and also a larger transverse flow for neutrons than 
for protons in heavy-ion collisions at higher energies, the 
neutron-proton differential flow is thus a measure of these 
two combined effec ts of the symmetry potentials on neu- 
trons and protons jl83l Il84| . 
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FIG. 25: (Color online) Rapidity distributions of the isospin 
asymmetry of free nucleons (upper panels), the difference of 
the average nucleon transverse flows (middle panels) and the 
neutron-proton differential transverse flow (lower panels) from 
132 Sn+ 124 Sn reaction at the incident beam energies of 400, 
800 MeV/nucleon and b = 5 fm with t wo sy mmetry energies 
of x = and x = -1. Taken from Ref. pM |. 

Shown in Fig. [25] are the rapidity distributions of the 
isospin asymmetry of free nucleons (upper panels), the 
difference of the average nucleon transverse flows (mid- 
dle panels) and the neutron-proton differential transverse 
flow (lower panels) from the 132 Sn+ 124 Sn reaction at in- 
cident beam energies of 400 and 800 MeV/nucleon and 
an impact parameter of b = 5 fm with the two sym- 
metry energies of x = and x = — 1. It is seen that 
the stiffer symmetry energy (x — —1) leads to a larger 
isospin asymmetry of free nucleons (stronger isospin frac- 
tionation) (upper panels) than the softer symmetry en- 
ergy (x = 0). As a result, the neutron-proton differential 
transverse flow from the stiff symmetry energy (x = —1) 
(bottom panels) is much larger than that from the soft 
symmetry energy (x = 0) compared to the difference be- 
tween the average neutron and proton transverse flows 
obtained from the two symmetry energies (middle pan- 
els). Also, the negative (positive) values of the neutron 
and proton average transverse flow F x — F p at forward 
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(backward) rapidities as a result of the Coulomb poten- 
tial effect on protons is reversed for the neutron-proton 
differential flow F^_ p (y) after taking into account the ef- 
fect due to isospin fractionation. 

From panels (c) and (f) of Fig. it is seen that the 
slope of the neutron-proton differential transverse flow 
around the mid-rapidity obtained from the same sym- 
metry energy is larger for the higher incident beam en- 
ergy. This is mainly because a denser nuclear matter 
is formed at higher incident beam energy, leading to a 
stronger symmetry potential and thus higher transverse 
momenta for neutrons compared to protons. Although 
the net effect of the symmetry potential on the neutron- 
proton differential transverse flow at 800 MeV/nucleon is 
not much larger than that at 400 MeV/nucleon, its mag- 
nitude is much larger and is thus easier to be measured 
experimentally. 
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FIG. 27: (Color online) Rapidity distribution of the double 
neutron-proton differential transverse flow in semi-central re- 
actions of Sn+Sn isotopes at incident beam energies of 400 
and 800 MeV/nucleon with two symmetry energies of x = 
and x = — 1. Taken from Ref. Il88ll. 



FIG. 26: (Color online) Same as panels (c) and (f) of Fig. I2"5l 
but for the reaction system of 112 Sn+ 112 Sn. Taken from Ref. 

EH. 

The systematic errors in the neutron-proton differen- 
tial flow can be reduced by studying its values in two 
similar reaction systems. Besides the 132 Sn+ 124 Sn re- 
action involving neutron-rich nuclei, one can consider 
another less neutron-rich reaction system 112 Sn+ 112 Sn. 
Fig. [221 shows the rapidity distribution of the neutron- 
proton differential transverse flow in the semi-central re- 
action of 112 Sn+ 112 Sn at the same incident beam en- 
ergies of 400 and 800 MeV/nucleon. In comparison 
with the 132 Sn+ 124 Sn reaction, the slope of the neutron- 
proton differential transverse flow around mid-rapidity 
and effects of the symmetry energy are much smaller 
due to the smaller isospin asymmetry in the reaction of 
112 Sn+ 112 Sn. 

The double neutron-proton differential flow is de- 
fined as the difference of the neutron-proton differential 
flows in the two reaction systems of 132 Sn+ 124 Sn and 
112 Sn+ 112 Sn. Fig. [2T1 shows the rapidity distribution of 
the double neutron-proton differential transverse flow in 



the semi-central reactions of Sn+Sn isotopes. At both 
incident beam energies of 400 and 800 MeV/nucleon, the 
double neutron-proton differential transverse flow around 
mid-rapidity is essentially zero for the soft symmetry en- 
ergy of x = 0. It displays, however, a finite slope with 
respect to the rapidity for the stiffer symmetry energy 
of a; = — 1. Moreover, the double neutron-proton dif- 
ferential transverse flow at the higher incident energy 
exhibits a stronger symmetry energy effect as expected. 
Since the double neutron-proton differential transverse 
flow retains about the same symmetry energy effect as 
the 132 Sn+ 124 Sn reaction, it is less s ensit ive to the sys- 
tematic uncertainties in experiments [158j |. 

Also, the Coulomb effect, which competes strongly 
with the symmetry potentials, is less important in the 
double neutron-proton differential transverse flow than 
in the neutron-proton differential transverse flow. This 
can be seen in Fig. [55] which shows the neutron-proton 
differential transverse flow (upper two panels) and the 
double neutron-proton differential transverse flow (low- 
est panel) in the semi-central reactions of Sn+Sn isotopes 
at the incident beam energy of 400 MeV/nucleon with the 
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FIG. 28: (Color online) Coulomb effects on the neutron- 
proton differential transverse flow (upper two panels) and 
the double neutron-proton differential transverse flow (low- 
est panel) in the semi-central reactions of Sn+Sn isotopes at 
the incident beam energy of 400 MeV/nu cleon with the sym- 
metry energy of x — 0. Taken from Ref. [I88t |. 



symmetry energy of x = for the two cases of with and 
without the Coulomb potential. From the upper two pan- 
els of Fig. [3H1 one sees that the Coulomb effect reduces 
the strength of the neutron-proton differential transverse 
flow as it makes more protons unbound and to have large 
transverse momenta in the reaction-plane. The Coulomb 
effect is, however, largely reduced in the double neutron- 
proton differential transverse flow shown in the bottom 
panel of Fig. El 



F. The n/p ratio of squeezed-out nucleons 

It is well known that in noncentral heavy-ion collisions 
nucleons in the participant region is squeezed out of the 
reaction plane as a result of the large density gradient 
in this direction and the absence of spectator nucleons 
to block their emissions. These nucleons can thus carry 
direct information about the high density phase of the 
reaction and have been widely used in probing the EOS of 
dense matter, particularly that of the symmetric nuclear 



matter, see, e.g., Refs. [| HH, HM GH, EM EH for a 

review. 
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FIG. 29: (Color online) Azimuthal distribution of midrapidity 
nucleons emitted in the reaction of 132 Sn+ 124 Sn at an incident 
beam energy of 400 MeV/nuc leon and an impact parameter 
of b = 5 fm. Taken from Ref. [Tol. 



The possibility of using the squeezed out nucleons to 
study the high density behavior of the nuclear symme- 
try ene rgy has been studied recently using the IBUU04 
model [l04| . Shown in Fig. [29] are the azimuthal dis- 
tributions of free nucleons in the midrapidity region 
(\(y/ybeam)c. m .\ < 0.5) of the reaction 132 Sn+ 124 Sn at 
an incident beam energy of 400 MeV/nucleon and an 
impact parameter of b = 5 fm, predicted by the IBUU 
model using the MDI interactions with x = and — 1. 
A preferential emission of nucleons perpendicular to the 
reaction plane is clearly observed for both neutrons and 
protons. The sensitivity to the nuclear symmetry energy 
is mainly seen in the squeezed-out neutrons, which expe- 
rience a stronger repulsive symmetry potential for a stiff 
[x = —1) than for a soft {x = 0) nuclear symmetry en- 
ergy in the neutron- rich matter, in addition to the strong 
nuclear isoscalar potential. For protons, their azimuthal 
distributions is less sensitive to the stiffness of the nu- 
clear symmetry energy. This is due to the additional 
repulsive Coulomb potential which works against the at- 
tractive symmetry potential experienced by protons in 
the neutron-rich matter. Although it is not easy to mea- 
sure neutrons in experiments, both the transverse flow 
and squeeze-out of neutrons together with other ch arged 
particles were measured at both the BEVALAC 
and SIS/GSI [HI EU Gil. These experiments and 
the assoc iated theoretical calculations, see, e.g., Refs. 
jl93l Il94| . have, however, all focused on extracting in- 
formation about the EOS of symmetric nuclear matter 
without paying much attention to the effects due to the 
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tides most likely come from the high density region in 
the early stage during heavy-ion collisions, they are thus 
more sensitive to the high density behavior of the sym- 
metry energy. Without the cut on the azimuthal angle, 
the n/p ratio of free nuclcons in the midrapidity region is 
much less sensitive to the symmetry energy in the whole 
range of transverse momentum, is as shown in the upper 
window. It is worth mentioning that the n/p ratio of free 
nucleons perpendicular to the beam direction in the CMS 
frame in 124 Sn+ 124 Sn reactions at 50 MeV/ nucleon was 
recently measured at the NSCL/MSU flQ5j | . This mea- 
surement was useful for studying the density dependence 
of the symmetry energy at sub-normal densities. To in- 
vestigate the symmetry energy at supra-normal densi- 
ties, similar measurements that allow the construction 
of the reaction plane using a TPC (Time-Projection- 
Chamber) and simultaneous detection of neutrons to- 
gether with cha rged particles at much higher energies are 
being planned [196| . The results reviewed here provide 
strong scientific motivations and support for such exper- 
imental efforts. Compared to other potential probes, the 
n/p ratio of squeezed-out nucleons is complementary but 
carries more direct information about the symmetry en- 
ergy at high densities. The sensitivity to the high density 
behavior of the nuclear symmetry energy observed in the 
n/p ratio of squeeze-out nucleons is probably the high- 
est found so far among all observables studied within the 
same transport model. 



FIG. 30: (Color online) Transverse momentum distribution 
of the ratio of midrapidity neutrons to protons emitted in the 
reaction of 132 Sn+ 124 Sn at the incident beam energy of 400 
MeV/nucleon and impact parameter of b — 5 fm, with (lower 
window) and without (upper window) an azimuthal angle cut 
of 80° < <t> < 100° and 260° < (j>< 280°, which corresponds to 
squeezed-out nucleons emitted in the di recti on perpendicular 
to the reaction plane. Taken from Ref. [I04j[ . 

To reduce the effect due to uncertainties associated 
with the EOS of symmetric nuclear matter, it is useful 
to consider the ratio of squeezed-out neutrons to protons, 
particularly its transverse momentum dependence. As 
shown in Refs. [J,[68j], the n/p ratio is determined mostly 
by the density dependence of the symmetry energy and 
almost not affected by the EOS of symmetric nuclear 
matter. In Fig. 1301 we show the transverse momentum 
dependence of the neutron/proton (n/p) ratio of midra- 
pidity nucleons emitted in the reaction of 132 Sn+ 124 Sn 
at the incident beam energy of 400 MeV/nucleon and 
impact parameter of b = 5 fm. For squeezed-out nucle- 
ons emitted in the direction perpendicular to the reaction 
plane, which are obtained by introducing an azimuthal 
angle cut of 80° < c/> < 100° and 260° < < 280°, the 
symmetry energy effect on the n/p ratio increases with 
the increasing transverse momentum p t as shown in the 
lower window. The effect can be as large as 40% at a 
transverse momentum of 1 GeV/c. Since high p t par- 



G. K°/K+ and E"/E+ ratios 

Since the proposal of Aichelin and Ko that the kaon 
yield in heavy ion collisions at energies that arc below 
the threshold for kaon production in a nucleon-nucleon 
collision in free space may be a sensit ive p robe of the EOS 
of nuclear matter at high densities |l97j | , a lot of works 
have been don e bo t h th e oretically and exper imentally on 
this problem [HE GH, GH, HOA, Sill IIqI. Since the 
kaon is an iso-doublet meson with the quark content of 
ds for K° and us for K + , the K°/K + ratio provides a 
potentially good probe of the nuclear symmetry energy 
as the n/p and 7r~ /tt + ratios, especially its high density 
behavior as kaons are produced mainly from the high 
density region during the early stage of the reaction and 
suffer negligible absorption effects. 

Using the UrQMD model (version 1.3), Li et al. have 
investigated the symmetry energy effects on the K°/K + 
ratio by studying K a and K + production from the central 
132 Sn+ 132 Sn collisions at a beam energy 1.5^4 GeV with 
two different forms of the symmetry energy, namel y, th e 
F15 and Fa3, and the results are shown in Fig. [31] [2031 ]. 
Since the beam energy is close to the kaon production 
threshold, which is about 1.58 GeV for nucleon-nucleon 
collisions in free space, the K°/K + ratio displays only 
a small symmetry energy effect. With decreasing beam 
energy, the symmetry energy effect becomes larger. For 
example, the kaon yields from the reaction 208 Pb+ 208 Pb 
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FIG. 31: (Color online) Top: Time evolution of the K° and 
K + abundances for central 132 Sn+ 132 Sn collisions at a beam 
energy 1.5 A GeV and with the symmetry potentials F15 
and Fa3. Bottom: the corresponding time evolution of the 



K°/K + ratios. Taken from Ref. [20: 



at E h = 0.8 A GeV and b = 7 ~ 9 fm, the K°/K+ 
ratio for the stiff F15 is about 1.25, whereas it is about 
1.4 for the soft Fa3. We note nuclear in- medium effects 
on kaon production are neglected in the UrQMD model 
simulations. 

Above results were confirmed by Ferini et al. using 
a relativistic hadronic transport model of Boltzmann- 
Uehling-Uhlenbeck type (RBUU) with different forms of 
the symmetry energies in ce ntra l (b = fm impact pa- 
rameter) Au+Au collisions 102]. Their results, shown 
in Fig. [321 indicate that at beam energies below and 
around the kinematical threshold of kaon production, the 
K° /K + inclusive yield ratio is more sensitive to the sym- 
metry energy than the n~ /tt + , and subthreshold kaon 
production thus could provide a promising tool to ex- 
tract information on the density dependence of the nu- 
clear symmetry energy. 

Experimentally, the FOPI collaboration has reported 
recently the results on K + and K° meson production in 
44 Ru + 44 Ru and 40 Zr + |gZr collisions at a beam kinetic 
energy of 1.528 A GeV , me asured with the FOPI detec- 
tor at GSI-Darmstadt [204|. The measured double ratio 
(K + /K ) Ru /(K+/K ) Zr is compared in Fig. [3j to the 




FIG. 32: (Color online) n~ /n + (upper) and K + /K° (lower) 
ratios as a function of the incident energy for central (b — fm 
impact parameter) Au+Au collisions with the RBUU model . 
In addition, for Eb eam = 1 AGeV, NLp results with a density 
dependent p-coupling (triangles) are also presented. The open 
symbols at 1.2 AGeV show the corresponding results for a 
132 Sn+ 124 Sn collision, more neutron rich. N ote t he different 
scale for the n~ /ty + ratios. Taken from Ref. |l02j |. 



predictions of a thermal model and the RBUU transport 
model using two different collision scenarios and under 
different assumptions on the stiffness of the symmetry 
energy. From Fig. 1331 one can see a good agreement with 
the thermal model prediction and the assumption of a 
soft symmetry energy for infinite nuclear matter while 
more realistic transport simulations of the collisions show 
a similar agreement with the data but exhibit a reduced 
sensitivity to the symmetry term. We note that in the 
present RBUU calculations, the isospin dependence of 
the K + - and if°-nucleon potentials in the asymmetric 
nuclear medium has been neglected. Recently, Mishra et 
al. studied the isospin dependent kaon and antikaon op- 
tical potentials in dense hadronic matter using a chiral 
SU(3) model and their results indicate that the density 
dependence of the isospin asymmetry is appreciable for 
the kaon and antikaon optical potentials. On the other 
hand, subthreshold kaon production in heavy-ion colli- 
sions depends on some de tailed implementations of the 
transport model [20 li |20 j |. Therefore, extracting useful 
information on the high density behavior of the nuclear 
symmetry energy from subthreshold kaon production in 
heavy-ion collisions induced by neutron-rich nuclei needs 
further studies from both experimental and theoretical 
sides. 

Besides the K°/K + ratio, the ratio has also 

been proposed as a probe of the high density behavior 
of the nuclear symmetry energy based on the UrQMD 
model (version 1.3) calculations (97[. Shown in Fig. 1511 is 
the time evolution of the 7r~ /ir + ratios (left-hand side) 
and the ratios (right-hand side) calculated with a 
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DATA THERM. NL NLp NLp5 

Data Vs Models 

FIG. 33: (Color online) Experimental ratio 

(K + /K°) Ru /(K + /K°) Zr (star) and theoretical predic- 
tions of the thermal model (cross) and the transport model 
with 3 different assumptions on the symmetry energy: NL 
(circles), NLp (squares) and NLp<5 (triangles). The INM and 
HIC calculations are represented by open and full symbols, 
respectively (see text for more details). The statistic and 
systematic errors are represented by ve rtical bars and 
brackets, respectively. Taken from Ref. [204 ]. 



stiff symmetry energy 1 and a soft symmetry energy 
F£= 3 for the reaction 132 Sn+ 132 Sn at E beam = 1.5A, 
2.5A, 3.5A GeV and b = 2 fm, and 112 Sn+ 112 Sn at 
Eb = 3.5A GeV and 6 = 2 fm. It is seen that the 
E~/E + ratio is sensitive to the density dependence of 
the symmetry energy for neutron-rich 132 Sn+ 132 Sn col- 
lisions, but insensitive to that for the nearly symmetric 
112 Sn+ 112 Sn collisions. For 132 Sn+ 132 Sn at E b = 1.5 A 
GeV, the E~/E + ratio calculated with the stiff symmetry 
energy (F^~ ) is higher than the one with the soft sym- 
metry energy (F^ =3 ). As the beam energy increases, the 
E~/E + ratio falls and the difference between the E _ /E + 
ratios calculated with and F£ =3 reduces strongly. 

As the beam energy increases further, at Eb = 3.5A GeV 
the E~/E + ratio falls further but the difference between 
the E~/E + ratios calculated with F^~ l and F% =z ap- 
pears again, the E~/E + ratio with soft symmetry energy 
now becoming higher than that with the stiff one. For 
pions, the results indicate that the ratio n~ /ir + at high 
energies (as in the case with Eb = 3. 5 A GeV) becomes 
insensitive to the symmetry energy. The difference be- 
tween the E _ /E + ratio and the n~ /t: + ratio can be un- 
derstood from the fact that, like nucleons, E^ 1 hyperons 
are under the influence of the mean field produced by the 
surrounding nucleons, as soon as they are produced. The 
symmetry potential of hyperons thus play an important 
dynamic role and results in a strong effect on the ratio 
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FIG. 34: The ratios vr" /tt+ (left) and S _ /E + (right) for the 
collisions 132 Sn + 132 Sn (E b = 1.5A, 2.5A, and 3.5 A GeV; 
6 = 2 fm) and 112 Sn + 112 Sn (E b = 3.5A GeV; 6 = 2 fm), 
calculated with the different symmetry potentials F^ 1 and 



of the negatively to positively charged E hyperons. 



VII. SUMMARY AND OUTLOOK 

Heavy-ion collisions induced by neutron-rich nuclei 
provides a unique opportunity for investigating the prop- 
erties of the isospin asymmetric nuclear matter, espe- 
cially the density dependence of the nuclear symmetry 
energy. To extract useful information from these colli- 
sions, transport models have been found to be extremely 
useful. Applications of these models have helped us 
understand not only the isospin dependence of the in- 
medium nuclear effective interactions but also that of the 
thermal, mechanical and transport properties of asym- 
metric nuclear matter. These information, particularly 
the density dependence of the nuclear symmetry energy, 
are very important for both nuclear physics and astro- 
physics. Significant progress has been made in recent 
years in determining the density dependence of the nu- 
clear symmetry energy. Based on transport model cal- 
culations, a number of sensitive probes of the symmetry 
energy have been identified. In particular, the momen- 
tum dependence in both the isoscalar and isovector parts 
of the nuclear potential was found to play an important 
role in extracting accurately the density dependence of 
the symmetry energy. From comparison of results from 
the transport model with recent experimental data on 
isospin diffusion from NSCL/MSU, a symmetry energy of 
E sym (p) « 31.6(/9/po) 7 with 7 = 0.69—1.05 at subnormal 
densities, which corresponds to the isospin and momen- 
tum dependent MDI interaction with x — and — 1, has 
been extracted. This conclusion is consistent with those 
extracted from studying other observables such as the 
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isoscaling data and the neutron-skin thickness in 208 Pb. 
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FIG. 35: (Color online) Density dependence of the nuclear 
symmetry energy using the MDI interaction with x — and 
x = — 1 and ot her many -body theories predi ction s (data are 
taken from [H, [HI Ho|])- Taken from Ref. flui ]. 

Although considerable progress has been made in de- 
termining the density dependence of the nuclear sym- 
metry energy at sub-normal densities, probing the high 
density behavior of the nuclear symmetry energy remains 
a major challenge. Our current knowledge on the density 
dependence of the nuclear symmetry energy can be seen 
from Fig. 1351 where pr edict ions from several typical the- 
oretical model [Tol . [l8L |205| are compared with the phe- 
nomenological constraints that have been obtained from 
model analysis of experimental data. The constraints la- 
beled x = and x = — 1 were extracted from studying 
the isospin diffusion in the reaction of 124 Sn + 112 Sn at 

cam — 50 AMeV withi n an isospin and m omentum de- 
pendent transport model (9llll07l . ll08lll09l |. For this par- 
ticular reaction, the maximum density reached is about 
1.2po- Moreover, it was shown that the neutron-skin 
thickness in 208 Pb calculated within the Hartree-Fock ap- 
proach using the same underlying Skyrme interactions as 
the ones labeled x = and x = — 1 is c o nsist ent with the 
available experimental data jlldl . Illll Ill2j . The sym- 
metry energy labeled as FSU-Gold was calculated within 
a RMF model using an accurately calibrated parame- 



ter set such that it reproduces both the giant monopole 
resonance in 90 Zr and 208 Pb as wel l as the isovector gi- 
ant dipole resonance of 208 Pb [ll3j . We note that the 
constraint obtained from the isoscaling analysis i s als o 
consistent with the FSU-Gold and the x = case |l!4| . 
At present, these results represent the best phenomeno- 
logical constraints on the nuclear symmetry energy at 
sub-normal densities. 

Although all predicted nuclear symmetry energies 
shown in Fig. [35] are close to the existing constraints 
at low densities, they diverge widely at supra-normal 
densities, including those from the MDI interaction with 
x = — 1 and that with x = as well as the FSU-Gold. 
Since there are currently no experimental constraints on 
the high density behavior of the nuclear symmetry en- 
ergy, more work is thus needed. In particular, experimen- 
tal data including neutrons from reactions with neutron- 
rich beams in a broad energy range will be useful for 
studying the behavior of the symmetry energy at high 
densities. We have reviewed recent theoretical progress 
in identifying the observables in heavy-ion collisions in- 
duced by high energy radioactive nuclei that are sensitive 
to the high density behavior of the nuclear symmetry en- 
ergy. A plethora of potentially sensitive probes have been 
found, and they include the ix~ ratio, isospin frac- 
tionation, n-p differential flow, double n/p and n~ /n + 
ratio, double n-p differential transverse flow as well as the 
K°/K + and £~/£ + ratios. Studying these observable in 
future experiments at high energy radioactive beam facil- 
ities is expected to lead to significant constraints on the 
behavior of the symmetry energy at supra-normal den- 
sities. Since transport models are essential in extracting 
these constraints from the experimental data, the contin- 
uous development of a practically implementable quan- 
tum transport theory for nuclear reactions induced by 
radioactive beams is important but poses a great chal- 
lenge. 
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